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ABSTRACT 

We present a novel class of models for Type la supernova time-evolving 
spectral energy distributions (SED) and absolute magnitudes: they are each 
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modeled as stochastic functions described by Gaussian processes. The values 
of the SED and absolute magnitudes are defined through well-defined regression 
prescriptions, so that data directly inform the models. As a proof of concept, we 
implement a model for synthetic photometry built from the spectrophotometric 
time series from the Nearby Supernova Factory. Absolute magnitudes at peak 
B brightness are calibrated to 0.13 mag in the g-band and to as low as 0.09 mag 
in the z = 0.25 blueshifted i-band, where the dispersion includes contributions 
from measurement uncertainties and peculiar velocities. The methodology can 
be applied to spectrophotometric time series of supernovae that span a range 
of redshifts to simultaneously standardize supernovae together with fitting cos- 
mological parameters. 

Subject headings: cosmology:distance scale - methods:data analysis - stars: supernovae: general 



1. Introduction 



Type I supernovae were found to be good sta ndard candles with a low dispersion in 
visible- wavelength peak luminosities (jKowall Il968l ). As supernovae were further subclas- 
sified, Type la supernov ae (SNe la) w ere found t o be even b etter calibrated candles by 
accounting for temporal (|Phiilipslll993l ) and color (|Tripplll998l ) information transmitted in 
their light curves. Modern distance determinators fit supernova measurements to a time- 
evolving spectral energy distribution (SED) model with free parameters that adjust the 
light-curve shapes and colors. An individual s upernova's peak ab solute magnitude is de- 
rived directly from the fit (e.g. Ijha et al.l l2007l: iMande l et alJboill, MLCS2k2, BayeSN) or 



from a secondary relati on based on the fit parameters (e.g. iGuv et al 



2007 



2008 



Conlev et al. 



Burns et al.ll201ll . SALT2, SIFTO, SNooPy). 



There are challenges in empirically modeling the family of SN la SEDs and relating 
supernova observables with absolute magnitude. Temporal coverage and sampling, and 
SN-frame photometric coverage vary from supernova to supernova; seldom are there data 
of different supernovae taken at common phase and band that can be directly compared. 
Data are fit to SED models to transform these heterogeneous observations into parameters 
that have common meaning across all supernovae. Constructing these models requires 
physical insight, diligent scrutiny of data, and/or extreme flexibility in order to ensure 
capturing the parameters that describe SN la heterogeneity. Once those parameters are 
identified, selecting the functional form to model their relationship with absolute magnitude 
is non-trivial: the behavior is known to be non-linear and there is the difficulty of choosing 
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a function that leaves no residual bias over the full parameter space while not overfitting 
the realized data set. We therefore seek an alternative to choosing the explicit functions 
for the time-evolving SED and the parameter-absolute magnitude relationship. 

This article introduces an alternative approach: the evolving SED for a supernova plus 
foreground dimming are modeled as a Gaussian process: a function whose value at each 
input is known to lie within a Normal distribution but whose precise value is unknow n and 
whose values at different inputs may be correlated (jRasmussen Williams! 120061 ) . The 
simple form of the Normal distribution leads to analytic expressions for the likelihood of 
the data for the Gaussian process model, and for the predicted mean and covariance of 
the fluxes at any wavelength and epoch. Modeled as a Gaussian process, the description 
of the underlying SED of a single supernova is influenced by each datum and thus has 
more degrees of freedom than either of the SEDs underlying MLCS2k2 and SALT2. Here 
supernova absolute magnitudes are modeled as a function of the flux values at different 
epochs and wavelengths regressed using the first Gaussian process. In one case a linear 
model is used. In a second Gaussian processes are again used to model supernova absolute 
magnitudes as a non-linear function. Previous applications o f Gaussian processes in the 
analysis of supernova data are given in iHolsclaw et al.l (|2010l ) ; IShafieloo et ajj (|2012l ) . 



Our approach contrasts with the usage of Gaussian processes in lMandel et al.l ( 201ll ). 
We model the covariance between single-band absolute magnitudes of different supernovae 
with a parameterized function that depends on differences between light-curve shapes and 
colors. The absolute magnitude of a supernova is regressed from the absolute magnitudes of 
other supernovae. The distance is then inferred from the observed and absolute magnitude 
in that band; indepe ndent informat i on po tentially available in the distances from other 
bands is not included. iMandel et al.l (|201ll ) focus on supernova distance and its covariance 
with a suite of light-curve parameters. There independent elements of the covariance matrix 
are fit in training. The distance of a supernova is then inferred directly from its light curve 
parameters and the trained covariance matrix. 

The approach presented in this article can be used in the analysis of existing supernova 
data and can infor m the design of futur e experiments. To date, the Nearby Supernova 
Factory (SNfactorv. lAldering et al.ll2002l ) has obtained and processed spectrophotometric 
time series for a sample of 132 supernovae within the linear Hubble flow. These data 
can serve as a template sample to search for indicators of diversity in luminosity, and 
contribute to the low-redshift anchor when combined with available high-redshift SN la 
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data and future photometric surveys such as the Dark Energ y SurvejQ, LSSl| or high- 
redshift spectrophotometric sur yeys as espoused by th e Interim Science Working Group of 



the Joint Dark Energy Mission (jBaltav fc Moosll201Cll ). 



In ^J2j we present in detail how Gaussian processes can determine SN la distances from 
broad-band photometry based on training with a low-redshift spectrophotometric sample. 
The framework for analyzing spectrophotometric data to calibrate absolute magnitudes 
and/or simultaneously fit for cosmological parameters is sketched in We finish with 
conclusions in 35l 



2. Case Illustration: Inferring Luminosity from Broad-Band Photometry 

The methodology to determine absolute magnitudes from broad-band light curves 
is applied to synthetic photometry generated from low-redshift SNfactory flux-calibrated 
spectroscopy (spectrophotometry). Two steps are involved. First, the photometry in all 
bands of all supernovae is modeled as a Gaussian process. Based on that model, interpo- 
lated multi-band light curves are constructed for each supernova. Then, the light-curve 
shapes and colors derived from the interpolations are treated as parameters that are re- 
lated to absolute magnitudes using a second Gaussian process. Based on that model, the 
absolute magnitude of a supernova with given light-curve shapes and colors is predicted. 

The motivation is to calibrate absolute magnitudes for supernovae at redshift z ob- 
served in the Dark Energy Survey (DES) griz photometric system. We work in the su- 
pernova restframe, treating observations as if taken in the appropriate blueshifted DES 
filter set. Focusing the analysis for specific redshifts is made possible for the first time 
by our extensive set of spectrophotometric data and eliminates K-correction uncertainties 



(jKim et al.lll996l ). In this article we consider four redshifts and thus four filter sets corre- 
sponding to the restframe coverage of z = 0.00,0.25,0.50,0.75 supernovae. The absolute 
magnitudes in each of the blueshifted griz bands at the epoch of -B-band peak brightness 
are calibrated. 

The analysis procedure for each redshift of interest is summarized as follows: 



1. Synthetic photometry of SNfactory data is generated. (£ j2.2 



J http : //www.darkenergysurvey . org/ 
^http : //www. lsst . org/lsst/ 
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2. The synthetic light curves from all supernovae are modeled as a Gaussian process. 
Best-fit parameters for the model are determined. ( §2.3.ip . 

3. SN light curves for a grid of epochs and bands and their associated covariance matrix 
are constructed by interpolating between a supernova's photometry using the regres- 
sion feature of Gaussian processes. Subtracting the interpolated absolute magnitude 
at i3-peak for a fiducial band gives light-curve shapes and colors. For computa- 
tional convenience the shapes and colors are re-expressed in terms of their principal 
components. ( §2.3. 2|) 

4. The absolute magnitude at B-peak as a function of the light-curve shapes and colors 
is modeled either as a second Gaussian process or a linear function and trained. 

(3233} 

5. The trained procedure is applied to predict absolute magnitudes for an independent 
validation set. This is done four times, running with completely independent vali- 
dation sets and corresponding complementary training sets. We report the average 
weighted rms (wrms) and standard deviation of the four-fold validation runs. ( §2.4p 



2.1. Computing 

Certain steps of Gaussian process regression involve the inversion or decomposition of 
matrices having dimension equal to the number of data points n. To handle large data sets 
(thousands of points and up), we implemented a framework of Gaussian process algorithms 
especially for use with high-performance parallel architectures. However, we note that there 
are alternative strategies for mitigating the C(n 3 ) matrix inversion cost, such as selective 
pruning of data and sparse linear algebra methods, but exploring them is beyond the 
scope of the present paper. Our i mplementation uses th e well-tested ScaLAPACK parallel 



distributed linear algebra library ([Blackford et al.lll997l ) 



The code base has the flexibility to operate on multiple platforms. Development was 
performed using one or two processors of a 2.53 GHz Intel Core 2 Duo MacBook Pro. 
The full analysis was run using 144 processors of Hopper, NERSC's Cray XE6. A single 
training/validation analysis took 5 wall-clock minutes to complete. 

High-performance computing facilitated code development and training: we were able 
to quickly cycle through tests using different Gaussian process models and optimize param- 
eters in likelihood fits. Once training of the model is completed, determining the absolute 
magnitudes of the validation set is not computationally intensive. 
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2.2. Data 



The input to our analysis is the portion of the spectrophotometric data set obtained 
by the S Nfactory between 2004 and 2009 with the SuperNova Integral Field Spectrograph 



(SNIFS, lLantz et al.l 120041 ) that had been fully processed as of early 2011. SNIFS is a 
fully integrated instrument optimized for automated observation of point sources on a 
structured background over the full ground-based optical window at moderate spectral 
resolution (R ~ 500). It consie 



Bacon et al. 


1995 


2000. 


2001 



, a multi- filter photometric 
channel to image the field in the vicinity of the IFS for atmospheric transmission monitor- 
ing simultaneous with spectroscopy, and an acquisition/guiding channel. The IFS possesses 
a fully-filled 6'.'4 x 6'.'4 spectroscopic field of view subdivided into a grid of 15 x 15 spatial 
elements, a dual-channel spectrograph covering 3200-5200 A and 5100-10000 A simultane- 
ously, and an internal calibration unit (continuum and arc lamps). SNIFS is mounted on 
the south bent Cassegrain port of the University of Hawaii 2.2 m telescope on Mauna Kea, 
and is operated remotely. Observations are reduced using the SN f actory 's dedicated data 



reduction pipeline, similar to that p resented in §4 of iBacon et al.1 (|200ll ). A discussion of 



the so ftware pipeline is presented in lAldering et al.l (|2006l ) and is u pdated in 



(|2010i ) . A detailed description of host-galaxy subtraction is given in 



Scalzo et al. 



Bongard et al.1 (|201lh . 



The instrument characteristics and observing strategy enable the construction of spec- 
trophotometric time-series of supernovae sampled every 2 to 3 days up until ~ 25 days after 
maximum, and then 5 to 7 days out to ~ 45 days after maximum. The median signal-to- 
noise within 2.5-days of peak brightness is 10.2 per 2.4 A bin. In our sample, 83% of the 
objects have z > 0.03 where peculiar velocities of order 300 km s _1 are expected to con- 
tribute less than 0.07 mag to the absolute magnitude dispersion (in fact 98% have z > 0.01). 
At the other extreme, 90% have z < 0.08 so are insensitive to present uncertainties in Qm 
or the dark energy equation of state, w. 



The spectral time-ser ies used are corrected for Milky Way dust extinction (jCardelli et al 



1989 



Schlegel et al.lll998l ) but no attempt is made to correct for the effects of circumstellar 
or host galaxy dust extinction. Re ference observations when the supernova has faded away 

Synthetic photometry light curves 



facilitate host-galaxy subtraction (IBongard et al 



2011 



are fit using SALT2 (|Guv et al.ll20o¥ to obtain the date of B-b&nd maximum in order to 
assign a phase to each spectrum. For our analysis, each spectral time series is blue-shifted 
to rest-frame and the fluxes normalized to represent that observed at a common distance 
for all supernovae assuming a flat universe with £Im = 0.28, h = 0.7 — thus, ratios in flux 
between supernovae can be interpreted as ratios in luminosity or equivalently differences 



- 7- 



in absolute magnitude. 

Synthetic photometry and uncertainties are generated for each of the blueshifted DES 
filter sets when there is sufficient wavelength coverage, including the few cases when data 
from only one arm of the SNIFS spectrograph is available. Spectra from the two arms gen- 
erally provide coverage of DES gri for z = 0.00, griz for z = 0.25, and riz for z = 0.50, 0.75; 
these redshifts provide different restframe wavelength ranges with which to calibrate super- 
nova absolute magnitudes. Photon fluxes are calculated and are converted to magnitudes, 
as we find that the the dominant sources of Gaussian dispersion are multiplicative (e.g. 
flux calibration) rather than additive (e.g. detector dark current). As the magnitudes re- 
fer to supernovae normalized to be at a common distance we refer to them as absolute 
magnitudes; the zeropoint of the magnitude system is irrelevant for the purposes of this 
article. 



Application of SALT2 on Johnson-Cousins synthetic photometry of the full input data 
set calibrates absolute magnitudes to a weighted rms of 0.147 mag, similar to the 0.16 mag 



found in earlier data productions (jBailey et al.l 12003 : IChotard et al.1 12011 



2.3. Modeling Supernova Absolute Magnitudes 

2.3.1. Modeling Supernova Light Curves 

We begin by training a parameterized model for the multi-band light curves that make 
up our data set. For a single blueshifted filter set, the multi-band light curves underlying 
supernovae are modeled as arising from a Gaussian process, 

™>(t,X) ~ GP (m(t, \]m ),k m (t, \,t' , \'-,l km ,a km )) , (1) 

where mrtX) 1S the photometric magnitude at epoch t relative to i?-band peak in the filter 
indexed by A . (The definition of a Gaussian process and associated equations used in the 
analysis of this article are given in $A]) Subscripted (t, A) represent points whereas those 
in unsubscripted parentheses are function variables. This notation means that for a set of 
coordinates {(t,\)}, the magnitude values are drawn from a Normal distribution 

m {(t,x)} ~N({m(t,\;m )},K), (2) 
where the elements of the model covariance matrix are 



K{t,\),(t',\') — k m (t, A,t , A';/ fcm ,o-fc m )- 



(3) 
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Measurements of a realized supernova are then also described by a Gaussian process, 

m (t,A) ~ GP {m(t,X;fao),k m (t, X,t', X';l km ,cr km ) + n m (t, A, t', A'; <r„ m )) . (4) 



Given a set of measurements m* t A ^ , the likelihood that the data arises from a Gaussian 
process can be expressed analytically, as well as the mean and covariance of the probability 
distribution function of a finite set of underlying points mr^x) ( see j Sfl) - 

For our model, we choose as mea n function m(t, A; mp) the synthetic photometry of 
band A from the updated templates of lHsiao et al.l (|2007l ) offset by zeropoints mo for each 
band. The mean function could include a parameter for the date of peak B magnitude but 
the spectral time series we work with has been preprocessed and is already expressed in 
terms of days relative to peak determined by SALT2. A zero mean function would lead to 
expectations of zero flux in temporal regions lacking data; the Hsiao mean is used to avoid 
abrupt deviations toward zero where there are temporal gaps in the light curves, whose 
locations differ from supernova to supernova. 

The kernel represents the covariance of the underlying flux: we use the square expo- 
nential kernel 



k m (t,X,t',X';l km ,a km ) = al m {X)5\x> exp 



t-t' 



(5) 



where each band has its own normalization factor a km (A) but shares a common time-scale 
l km . The lk m does not represent the time scale of the light curves themselves, but rather 
those of deviations from the mean model. The kernel does not correlate the light curves 
in different bands. Although supernova flux across bands is correlated, for this stage of 
the analysis we preferred imposing no prior assumptions on color. We opt for the square 
exponential kernel as it is frequently used in the Gaussian process literature; model testing 
comparing the likelihoods of other kernels is deferred for future study. 

Measurements of m include a nugget term that contains all other sources of variance 
in the photometry. Our nugget model is 



i(t, A, t , A'; <r nn 



[a 



l m (X) + var(m (tjA) )) 6 t;t ,5 X) \', 



(6) 



where within the scope of our analysis there are no two measurements with the same (i, A) 
so the nugget applies per measurement. Each band has its own intrinsic dispersion and 
the SNfactory photon noise has variance var(m(f^)). The measurement covariance is very 
small given the small overlap of the griz bands, and so is ignored. 
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The kernel inc ludes an additional term that accounts for radial peculiar velocities 



(IBonvin et al J 12003 ): 

2 nec ^ (1 + Z a ) 

P = P In (10) X H \ (7) 

where the peculiar velocity dispersion, crf ec , can be either a fixed or fit parameter, and 
z a is the actual redshift of the supernova; in this article we fix <rf ec = 300 km s _1 . Here 
the a and j3 are supernova indices, so the 5 a p reflects the common realization of peculiar 
velocity for a single supernova. At low redshifts the conformal distance can be conveniently 
approximated as % ~ z — 0.75J1M-Z 2 for a flat ACDM cosmology. 

For each of four training/validation runs, 75% of the SNfactory supernovae are desig- 
nated as the training set. Photometry with S/N < 50 are culled as often being associated 
with poor extractions and a supernova must have at least 16 light-curve points total over all 
bands to be included. Cuts are applied after the division into training and validation sets, 
so each of the analyses do not have exactly the same sample sizes: after the signal-to- noise 
and number-of-point cuts there remains 80-90 training supernovae per run. 

All training-set supernovae are used to fit the free parameters (and hyperparameters) 
that give the maximum likelihood of the Gaussian process model. The hyperparameters 
h m > &k m > and a- nm are common to all SNe la and are fit simultaneously to the full ensemble 
of data, whereas a distinct mo is fit for each supernova. The best-fit hyperparameters 
for Zfc m , o"fc m , and cr nm for one of the 4-fold training runs are shown in Table [TJ these 
values do not differ significantly between the different realizations. The strength of the 
model covariance between points in a single light curve range from 0.08-0.14 mag with 
correlation-length scales of around 6 days. The intrinsic dispersion is 0.05-0.07 mag except 
in the restframe UV where the dispersion is 0.14 mag. 



System z 










°k m (z) 


o"n m (g) 




0"n m (*) 




0.00 


5.63 


0.09 


0.09 


0.12 




0.05 


0.06 


0.07 




0.25 


5.87 


0.13 


0.08 


0.09 


0.13 


0.14 


0.05 


0.06 


0.06 


0.50 


6.18 




0.11 


0.08 


0.10 




0.06 


0.05 


0.05 


0.75 


6.56 




0.10 


0.11 


0.11 




0.14 


0.06 


0.05 



Table 1: Best-fit hyperparameters of the Gaussian process that models the supernova multi- 
band light curves, given in Equation [TJ Given are the results of one of the four training 
sets for each of the blueshifted DES photometric systems. All values have magnitude units 
except lk m , given in supernova- frame days. 
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The regressed g-band date of maximum has a distribution with mean and standard 
deviation of 0.058 ± 0.76 days. The temporal alignment of our model is consistent with 
that of SALT2. 

To summarize this first step in the analysis, supernova light curves are described by 
a parameterized model for the covariance in their residuals about a parameterized mean 
function. Light-curve photometry of a supernova training set are used to optimize the 
parameters that best describe their scatter. 

2.3.2. Supernova Parameterization 

We now turn to the supernova parameters used to describe the light curve of a single 
event, in analogy to the light-curve shape and color parameters used in classical light- 
curve fitters. The parameters adopted in this article are the (linearly transformed) values 
of interpolated light curves. 

Since we are ultimately interested in determining the magnitude at peak brightness, 
both training and validation sets are pared down to include only those supernovae having 
a first measurement at least two days before B peak to avoid extrapolating the value of 
M(^ o to ). This leaves 49-57 training supernovae. 

Observed supernovae have disparate temporal sampling; those data must be trans- 
formed into a set of parameters that can be directly compared between supernovae. The 
Gaussian process model of the light curve is used to predict the fluxes on a common set 
of epochs, daily from —10 to 35 rest-frame days relative to B peak, in the blueshifted 
DES bands that are spanned by SNfactory data. (The choice of daily sampling of the 
light curves was arbitrary and probably finer than necessary.) The predictions (Eqn. IA5h 
have covariance (Eqn. IA6h so fifty light-curve grids are realized per supernova to include 
uncertainties in the training and propagate uncertainties in the inferences. Predictions and 
realizations of regressed light curves from sample supernova photometric measurements are 
shown in Figures [H [21 and 

We want to calibrate the supernova absolute magnitude in some band Ao at the phase of 
peak B brightness t$: this is referred to as the "true" absolute magnitude and is estimated 
by the M^o to \ predicted by our model. This information is removed from the predicted 
light curves by subtracting M^x 0: t ) from all other magnitudes. What remains is a light- 
curve shape in the selected band Ao and colors for the other bands, giving 137 (46 x 3 — 1 
for the case of three bands) or 183 (46 x 4 — 1 for four bands) potential parameters to 




Fig. 1. — Synthetic photometry, predicted mean (solid), and ten realized (dashed) griz 
light curves for a random validation supernova (not used in the training) in the run with 
z = 0.25 filter set. 
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Fig. 2. — Light curves for SN #19 in the z = filter set with an arbitrary normalization 
offset. The predicted peak magnitudes for this object are off by ~ 0.6 mag. 
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Fig. 3. — Light curves in the z = filter set for SN #13, whose predicted peak magnitudes 
are off by ~ 0.3 mag. 
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describe a supernova realization. 

For practical reasons, we re-express the light-curve shapes and colors in terms of 
their principal components ordered by the dimensions that carry the leading population 
variance by applying a principal component analysis (PCA) on all 50 realizations of all 
training supernovae. The coordinates of the leading eigenvectors that account for the 
95% of the variance are used as the supernova parameters, and are denoted as x. This 
reduces the dimensionality of the problem down to ~ 15 in anticipation that the truncated 
components carry little real signal that correlates with absolute magnitude. As a result of 
this transformation each parameter no longer represents a specific color at a specific phase, 
but rather the contribution to the temporal evolution of all colors from a single PCA 
eigenvector. Note that due to the correlated noise in the regressed light curves, PCA is not 
encapsulating pure supernova diversity efficiently in a minimal number of components. 

As a representative example of one of the training sets, we show results from one of 
the runs calibrating the i-band of the z = 0.25 filter set. The first several PCA eigenvectors 
describing light-curve shapes and colors (after subtraction of the mean) are shown in Fig- 
ure [H These eigenvectors capture 54%, 15%, 9%, and 5% of the population variance. The 
distribution of PCA coefficient values for the training set is shown in Figure [5j For each 
component i, the coefficients have standard deviation a Xi . In solid (red) are the coefficients 
for the mean light curves and colors for each training-set supernova. Surrounding each of 
those points is a cloud of gray points of the coefficients from the fifty realizations; the size 
of the cloud is reflective of the dispersion in the Gaussian process model as well as the 
data quality. The ~ 50 training-set supernovae do not fill a continuous distribution in this 
space, there are clear outliers due to realizations of a single supernova. The distinct blobs 
containing multiple supernovae seen in the x(0)-x(2) plot are not significant as they are 
not as prominent in the three other training runs. 

Mean and perturbed light curves from 1- and 2-a Xj changes in the x(j) parameter 
are shown for j = 0, 1, 2 in Figures [6l [8] and [9l The j = eigenvector (Figure [6|) behaves 
as if it were due to dust: the correction is (mostly) phase-independent and monotonically 
decreases from blue to redder bands. The j = 1 eigenvector (Figure [7|) generates light 
curves that vary in width, with supernovae with relatively skinny blueshifted-g (observer- 
frame central wavelength of 380 nm) light curves having relatively broader light curves in 
the redder bands. The g— X colors varies but the other colors do not. The j = 2 eigenvector 
(Figure [8]) has a particularly strong influence in the shape and colors in the near-UV band. 
It has the property that the curves intersect at later phases; the effect is pronounced in 
the blueshifted i band where the phase and height of the secondary maximum vary. Light- 
curve shapes and all color combinations also are influenced by this component. The j = 3 



Fig. 4. — First four PCA eigenvectors of the z = 0.25 filter set plotted as a function of 
corresponding phase relative to B peak split into g (top left), r (top right), i (bottom left), 
and z (bottom right), from a run calibrating the i flux (hence the gap in the bottom left 
panel). These eigenvectors capture 54%, 15%, 9%, and 5% of the population variance. 



- 16 - 




-1 



-4 -2 2 4 6 -3-2-1 1 2 3 4 -2 -1 1 2 
x(0) x(1) x(2) 

Fig. 5. — Distribution of the first four PC A coefficients for one of the runs calibrating the i- 
band of the z = 0.25 filter set. The predicted mean of each training-set supernova is shown 
in red with the 50 associated realizations shown in black. The ellipses have semi-major 
and semi-minor axes with the standard deviations of the distribution a x . . 
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eigenvector (Figure [9]) is similar to j = 2 and has a particularly exaggerated influence on 
the shape of the near-UV band. 




Fig. 6. — Mean and perturbed light curves from 1- and 2-a XQ changes in the x(0) parameter 
for the z = 0.25 filter set from a run calibrating the i flux, plotted as a function of 
corresponding phase relative to B peak. The four bands are plotted: g (top left), r (top 
right), i (bottom left), and z (bottom right). 

The PCA coefficients are correlated with the SALT2 color and x\ parameters. Fig- 
ure [10] contains scatter plots of the first four principal components with the SALT2 color 
and x\ parameters for one of the runs for the z = filter set: the first PCA component is 
strongly correlated with color and the second and fourth are correlated with x\. 

To summarize this second step of the analysis, the mean multi-band flux and covariance 
of the underlying light curve of an individual supernova are calculated for a set of epochs 
and filters using the sparsely measured photometry combined with the trained light-curve 
Gaussian process model. Each realization of the light curves is decomposed into one "true" 
absolute magnitude and light-curve shapes and color curves. Anticipating that absolute- 
magnitude variations are encoded in SN la heterogeneity, for the light-curve shapes and 
colors we transform to a new coordinate-system that is defined with a PCA on the training 
set. 




Fig. 7. — As in Figure [6] except for perturbations in the x(l) parameter. 



2.3.3. Parameterized Model for Supernova Absolute Magnitudes 

The light-curve and shape coordinates in the PCA basis are correlated with true 
absolute magnitude (modulo a constant offset for the arbitrary distance at which the SNe 
are placed, see ^2.2j) . as seen in Figure [TTJ 

In this section we relate light-curve and color parameters to absolute magnitudes with 
a linear and a Gaussian process model. The linear model makes an explicit choice for 
the relationship between light-curve parameters and absolute magnitude that can be well- 
constrained by our small dataset. The Gaussian process methodology lets the data drive 
the model. While having the advantage of not requiring an explicit model, results can be 
sensitive to statistical fluctuations in the training set or large interpol ations in sparsely 



populated regions. An analogous contrast can be made between SALT (jGuy et al.ll2005l ) 



who used parameterized functions, and MLCS (jRiess et al.lll996l ) who used data averages 
to construct light curve models. 
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Fig. 8. — As in Figure [6] except for perturbations in the x(2) parameter. 

Linear Model 

The absolute magnitude at phase to an d band Ao is a linear function of the light-curve 
parameters M( to! ^ )(x;p^) = Mq + Y^]=o x ~ 1 Pj x U)i where Mq and pj are fit parameters. 
Most of the over 100 supernova parameters capture little of the population variance: as is 
justified in £ 12.41 n max = 4 is adopted in this article. 

Figure [T2l shows the linear-model solution for the z = 0.00-filter set calibrating the g 
band, which closely corresponds to the bands for which SALT2 is directly trained. As in 
SALT2, there are clear correlations between the first two parameters, which in turn are 
correlated with color and x\ respectively (as shown in Figure [T0|) . There is no correlation 
for the third parameter, but a correlation with the fourth parameter is apparent. 

Gaussian Process Model 



While previous analyses have modeled linear and quadratic correlations between light-curve 
parameters and absolute magnitude, it is well known t hat such simple relat ionships are not 
representative of the full range of light-curve shapes (jHamuv et al.lll996l ). Here we infer 




Fig. 9. — As in Figure [6] except for perturbations in the x(3) parameter. 
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Fig. 10. — Scatter plots of the first four PCA coefficients, x(0), x(l), x(2), and x(3) with 
the SALT2 color and x\ (light-curve shape) parameters for one of the runs for the z = 0, 
filter set. The correlation coefficients of color with the PCA coefficients are {—0.91. — 
0.33, 0.02, -0.03} and of xi with the PCA coefficients {-0.18.0.80, 0.27, -0.47}. 
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Fig. 11. — Training-set absolute magnitudes (with constant offset) as a function of the first 
four light-curve coordinates in the PCA basis for one run of the z = 0.25-filters calibrating 
the i band. Each point represents one supernova realization. 
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Fig. 12. — In red the 1-er range of the inferred absolute magnitude found for the linear fit 
and in black the 1-a bands from the Gaussian process model as a function of parameter 
value while holding the coordinates for all other dimensions at zero, for one run of the 
z = 0.00-filters calibrating the g band. Note the different ordinate ranges of the plots. 
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the absolute magnitude through regression in the high-dimensional PCA parameter space. 

The fiducial magnitude at phase to and band Ao, M( tQ ^ ), is modeled as a Gaussian 
process dependent on light-curve shapes and colors expressed by the parameters x, 

M (toM) ~ GP ( M ( x ; M o), k M (*, x'; l kM , a kM ) + n M (x, x'; cr„ M )) . (8) 

The mean function M(x;p^) = Mq + ^""q 1-1 Pjx(j) is linear and has free parameters 
Mo and pj. The kernel function is again a square exponential 

k M (x, x'; l kM ,a kM ) = a\ M exp 

where n is the number of elements in x, i.e. the number of eigenvectors capturing 95% of the 
population variance. There is an independent per-point nugget, um(x-, x'; cr„ M ) = cr^ lM 5 x y. 

Training-set supernovae are used to find the best-fit values for Mo, l kM , Pj, &k M -> and 
cr nM . The values of cr ku , the strength of magnitude correlations, are found to be around 
0.18 mag. The residual absolute magnitude dispersion after training is given by the cr nM 's, 
which range from 0.04-0.07 mag except for the z = 0.75 bands that have a range from 
0.07-0.11 mag. 

The hyperparameters for the length scale of the correlation \ kM alone do not provide 
much physical insight. A more interesting statistic is the ratio of the correlation scale and 
a typical value of the parameter in the training set; the metric represents a typical scale 
for each parameter in the exponent in Eqn. [9l which in turn gives the relative contribution 
to the correlation strength. Parameters with low values in Fig. [13] contribute more to the 
correlation. Figure [TBI shows the best-fit l kM (j ') /cr Xj , where a x . is the standard deviation 
of the j'th coefficient of all realizations in the training set of a representative run. 

The relationship between each parameter and absolute magnitude is illustrated for a 
slice in parameter space for the same run shown for the linear model; Figure [J_2] shows 
the most-likely absolute magnitude as a function of parameter value while holding the 
coordinates for all other dimensions at zero. In this example the coordinate range spanned 
by the supernova set decreases with the dimension order: —6 < x(0) < 4, —2 < x(l) < 2, 
— 1.8 < x{2) < 1.8, and —1 < x(3) < 1. The correlation between parameters and absolute 
magnitude is non-linear. 

The Gaussian process and linear solutions for the first two parameters are in general 
agreement, consistent with the SALT2 model. X(2) does not have a monotonic trend in 
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Fig. 13. — Characteristic correlation scale in the kernel for each parameter, {i)/(T Xi , 
where o~ x . is the standard deviation of the i'th parameter over the training set. 

the region and extends a small range of influence on luminosity. The fourth parameter has 
agreement in X(3) within the range where there is training-set data. 

To summarize this step of the analysis, supernova absolute magnitudes are taken to 
be a function of the light-curve shapes and colors. In one analysis a linear fit is performed. 
In the other the functional form is not explicit, instead its stochastic behavior is modeled; 
training-set supernovae provide several measurements of this function and we use these to 
train a Gaussian process description of their correlation. The trained Gaussian process 
provides a prescription for interpolating an "inferred 1 absolute magnitude from the light- 
curve shape and color. It should be emphasized that the training procedure is not optimized 
to minimize the dispersion in the residuals of absolute magnitudes but rather to maximize 
the likelihood of the data for the Gaussian process model. 



2.4. Validation 



The validation-set supernovae are used to test how well true absolute magnitudes can 
be predicted from the inferred absolute magnitudes derived from the trained model in the 
procedure given in §2.3.11 We perform cross-validation by dividing the supernovae into 
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distinct training and validation sets with 25% assigned to the latter: the results are not 
expected to be as good as if all supernovae are used for training, but in this approach there 
is no question of the independence of the magnitude dispersion of the validation set from 
the data used in training. We measure the stability of our results by repeating the same 
analysis on validation sets containing different supernovae; as 25% of the objects enter 
the validation set the analysis is run four times with independent validation but partially 
intersecting training sets. Over four runs each supernova enters a validation set once and 
only once. Comparing the four almost-independent cross-validations diagnoses whether the 
number and fraction of training supernovae are a fair representation of the full population. 



2.4-1- Pruned Supernova Sample 

The SNfactory data do not have perfect temporal sampling or signal-to-noise. As a 
selection criterion, we study those validation supernovae with synthetic photometry in at 
least eight distinct phases in each of the bands. This leaves 13-23 validation supernovae 
per run except in the z = 0.75 filter set, which has 7-11 SNe per run. Measurement uncer- 
tainties contribute to uncertainties in the interpolated absolute magnitude and light curves. 
To illustrate, the distribution of light curves realized from the Gaussian process predicted 
mean and covariance shown in Figures [IJ [2j and [3] indicate the size of the uncertainties. 
Figure [15] contains a histogram of the standard deviations of the measured true absolute 
magnitude at -B-peak for the i-band of the z = 0.25 filter set, here the absolute magnitude 
uncertainties per supernova range from 0.02-0.06 mag. 

As a representative example of the light curve regression, we present results for one 
random supernova from one of the runs with z = 0.25 filter set. Figure [T] shows the 
synthetic photometry, the predicted mean of the underlying light curves in solid, and ten 
random realizations in dashed lines. 

We identify two odd supernovae, SN #13 and SN #19, that are culled from the sample. 
SN #19 data and predicted light curves for the z = filter set are shown in Figure [2j This 
supernova does not have the wavelength coverage to make the z = 0.25 sample and its 
PCA coefficients fail the distance cut for z = 0.75. It is an extremely slow decliner with a 
<7-band Amis = 0.7 derived from the Gaussian process prediction. The PCA coefficients 
from the z = case are shown among those of the training set in Figure [T3J it is clear 
that SN #19 does not have a corresponding analog in the training set. With z = 0.0055, 
this supernova has the lowest CMB-frame redshift of the full sample of which > 96% has 
z > 0.015; a peculiar velocity of 300 km s _1 would contribute a 0.40 mag residual. SN #13 
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data and predicted light curves for the z = 0.25 filter set are shown in Figure [3l The 
discrepancy is accounted for by one spectrum near peak that is fainter than expected; 
indeed this discrepancy disappears in subsequent versions of the spectral reduction. 




Fig. 14. — Distribution of the first four PCA coefficients for one of the z = 0.00 filter- 
set runs used in the training of the analysis of SN #19. The symbols are the same as 
in Figure [5l In addition, blue stars show the PCA coefficients of the predicted mean of 
SN #19. 

Validation supernovae may not have an analog in the training set. We approximate 
the PCA coefficient distribution as Gaussian, described by the standard deviations a Xi of 
each dimension. A validation supernova is included in the analysis if it lies within the core 
of the distribution 

Ef— Y <X 2 (p,np C A), (10) 

where x 2 {Pi n ) ls the value where the cumulative x 2 -distribution has p- value 0.05 and npcA 
is the number of principal components that capture 95% of the variance. This cut criterion 
was made after the distributions were constructed and is therefore not blind. However, 
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Fig. 15. — Histograms of uncertainties in the Gaussian process prediction variance including 
measurement error in the true (unfilled) and inferred (hatched) absolute magnitudes per 
supernova, for a run calibrating the i band of the z = 0.25 filter set. 

0.05 is selected as a nice round number that should capture 95% of the population and 
ensuing results are insensitive to departures from this value. This cut typically removes 
0-2 objects from the sample. 



2.4-2. True and Inferred Magnitudes and their Differences 

The covariance in the regression is accounted for in our analysis by using fifty light- 
curve realizations of each supernovae. Each of the fifty multi-band light-curve realizations 
of each supernova is temporally sampled from —10 to 35 days using the best-fit Gaussian 
process model of the training set. The light curves are normalized by the true Mu Qi x ) and 
expressed in the coordinate system defined by the PCA eigenvectors. The PCA coefficients 
are used to infer Mu Q ^ \- The standard deviation of its fifty realizations provides the 
measurement uncertainty of M^ t(hXo ) of a supernova, this value depends on the quality of the 
measured light curve and ranges. The means and standard deviations of the 50 realizations 
of true and inferred absolute magnitudes from the validation procedure (modulo a constant 
offset) are given in Tables \M5\ El and [7] for the Gaussian-process-model. (In the tables, 
omitted supernovae were used for light-curve training but did not pass the selection criteria 
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for the peak-magnitude analysis. In italics are supernovae that have data before two days 
pre-maximum with light curves that have at least eight points. In bold are supernovae used 
for validation that have data before two-days pre-maximum and at least eight points in 
all light curves.) Figure [T5l shows the histogram of the standard deviations of the inferred 
absolute magnitude at B-peak for the i-band of the z = 0.25 filter set: the inferred absolute 
magnitude measurement uncertainties per supernova mostly fall within the range 0.04 - 
0.10 mag. 

The accuracy to which true absolute magnitudes are predicted by absolute magnitudes 
inferred from SNfactory synthetic light curves is the root-mean-square of the difference be- 
tween the true and inferred M^ x ) over all realizations of all supernovae. The weighted 
rms (wrms) of the differences of each validation set, and then the mean and standard 
deviation of those weighted rms' over the four training/validations (wrms and a(wrms) 
respectively) are calculated. The wrms weights the squared residual by the inverse mea- 
surement variance y isualized in Figure [15] (the expression for the weighted rms is given in 



Blondin et al.ll201ll ). The mean shows how well the algorithm standardizes the distance 
indicator over all supernovae. The standard deviation of the wrms of the four runs indi- 
cates the sensitivity to the sample selection. The weighted rms' have contributions from 
measurement errors, peculiar velocities, and the intrinsic ability that the methodology has 
to calibrate supernovae, Oi n %. Assuming a peculiar velocity dispersion of 300 km s -1 , for 
each run we report the mean and standard deviations of the fit value the intrinsic dispersion 

The wrms calculated above represents an estimate of the absolute-magnitude uncer- 
tainty when using 75% of the data training, such that none of the validation dispersions 
are calculated with supernovae used in the training. Our analysis lends to easy calculation 



of the wrms from K-fold cross analysis (as done by lBlondin et al.ll201ll ) for K = 4, where 
the dispersion is calculated for all supernovae together despite the fact they are based on 
different training sets; these results are given in Tables [2] and El We provide two numbers 
for the K-fold analysis; the first using the subset after applying the sample cut of Eqn. [10] 
and the second using the full sample of objects. The former statistic is calculated antic- 
ipating that the supernova parameter space is not densely sampled and considering that 
our absolute magnitude model may not extrapolate well to objects outside the training 
regime. 

For the linear model, the number of light-curve parameters related to absolute magni- 
tude is set as n max = 4. This choice is made as follows: We calculate wrms for the linear 
model varying n max from 1 to 8 for the i-band calibration of the z = 0.25 filter set, getting 
0.097, 0.099, 0.104, 0.091, 0.092, 0.094, 0.093, and 0.095 mag. There is a step in the wrms 
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when four parameters are used. In the Gaussian process model, n max is set to the number 
of principal components that capture 95% of the variance. 

Calculated dispersions from the validation are shown in Table [2] for the linear fit, and 
in Table [3] for the Gaussian process model. The weighted mean of each validation set is 
< 0.01 mag. 



The statistics in lEfronl (|1983l ) are calculated for the case of the i-band from the z = 0.25 
filter set: The apparent error rate is determined from the residuals of the training-set super- 
novae as run through the same processing as the validation supernovae. The linear model 
gives 0.089 mag. The resulting 0.025 mag for the Gaussian process model is an overopti- 
mistic error estimate: the best-fit a nM = 0.060 mag nugget parameter (introduced after 
Eqn. [9]) gives an estimate of the absolute magnitude dispersion of the training supernovae. 
The bootstrap estimate is calculated from 50 realizations of training sets generated using 
random draws with replacement from the full set and the ~ 37% unselected supernovae 
forming the validation set; this error estimate is 0.107 (0.117) mag for the linear (Gaussian 
process) model. The .632 estimator is not rigorously motivated but empirically exhibits 
small biases in the estimator of the true prediction error in certain cases, it is a function 
of the apparent and bootstrap estimators for which we get 0.100 (0.083) mag for the linear 
(Gaussian process) model. 



3. Discussion 

Training and validation sets have different variance in their predicted magnitudes. 
The discrepancy between the two is not surprising given the sparsity of supernovae; it 
indicates that the current training set does not fairly represent the population of the 
current validation set. We anticipate that these values should converge as the number 
of training supernovae increases. The validation wrms may decrease with interpolation 
from more local values, while the training standard deviation may increase as singletons 
disappear. 

The subsample of ~ 50 training-set supernovae used to calibrate absolute magnitudes 
is too small to cover the full range of intrinsic and extrinsic supernova variability. This issue 
is highlighted in SN #19, where highly-discrepant absolute magnitudes are predicted when 
no analogue exists in the training set. We re mark that SN #19 (IAU name SN 20071e) 



does have unique properties noted elsewhere (jSimon et al.ll2009l ): it is one of only three 
supernovae that exhibit time- varying Na I D absorption features and the only one of those 
that is not highly reddened, suggesting circumstellar extinction. Nevertheless, with the 
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Table 2: Results from the linear model: For each filter system and calibrated band: <7o, the 
raw standard deviation of M( to Ao ) for the input sample; wrms and a(wrms), the mean and 
standard deviation over the four validation samples of the wrms difference between true and 
inferred magnitude; the mean intrinsic dispersion and its standard deviation over the four 
validation samples, the 4-fold cross-validation error estimate with (and without) rejection 
of objects with extreme parameter values; N$n is the number of validated supernovae. 
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Table 3: As in Table [2] for the Gaussian process model. 
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loose requirement of Eqn. [10] many other supernovae manage to be reasonably calibrated. 
With a larger sample size, it will be interesting to see if large outliers persist despite having 
training-set analogues in regions of parameter space where SNe la are not well standardized. 
The trade between criteria for deciding if a supernova has training-set analogues versus 
the numbers of supernovae that meet that criteria should be studied in actual cosmological 
analyses. 

Whether the linear or Gaussian process model results in lower dispersion depends on 
the filter set being calibrated and the specific statistic of interest. Figure [T2l shows examples 
of how either model can go wrong. The wiggles in the X(l)- dependence of absolute 
magnitude suggest that the Gaussian process regression is sensitive to fluctuations in the 
small training set. The Gaussian process regression for X{2) is clearly not well described 
by the linear model. While both linear and Gaussian process solutions for X(3) agree 
in the region covered by the training sample, beyond they extrapolate out with different 
slopes. For the data used, neither model for absolute magnitude universally outperforms 
the other. 

We get dispersions of 0.09-0.14 mag (linear fit) and 0.09-0.11 mag (Gaussian process 
model) using the four bands of the z = 0.25 filter set. Filter sets that use only three bands 
exhibit higher dispersion, although only modestly so. Although the input magnitude dis- 
persion strongly depends on the band of the magnitude being calibrated, this dependency 
mostly disappears after calibration in the Gaussian process fit. Some of the dispersion is 
due to limitations in the data themselves. There are light curves with suspected misesti- 
mated fluxes or uncertainties, SN #13 being a prime example (Figure [3j). Sometimes there 
are significant temporal interpolations being made due to coarse light curve sampling, par- 
ticularly over the secondary maximum in the redder bands as seen in Figures [IJ [21 and 
A further decrease in the dispersion may be possible as more supernova spectrophotometry 
with higher signal-to-noise and denser temporal sampling becomes available. 

In an attempt to avoid an explicit description of the supernova light-curve model, we 
use the relatively simple normalized Hsiao templates for the mean function. Certainly, 
Gaussian processes can be used to describe residuals about the SALT2 model whose shape 
and color parameters (including date of maximum) are included as variables in the absolute 
magnitude model. Inclusion of a dust-absorption model in the mean function should reduce 
the number of training-set supernovae needed to densely sample the full range of dust 
absorption. 

The new expanded set of light-curve parameters can standardize SNe la to a level that 
appears competitive with other fitters. A direct comparison with the quoted dispersions 
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of 0.12 mag for SALT2 and 0.14 mag of MLCS2k2 is difficult because there is no sharp 
delineation between training and validation supernovae reported for these fitters. A SALT2 
analysis on synthetic photometry from this same dataset yields a 0.147 mag residual on the 
Hubble diagram. An altern ative reference is th e 0.18-0.22 mag scatter found using these 
fitters on the CfA3 sample (IHicken et al.ll2009l ). Assuming the dust parameter Ry = 3.1, 
Mandel et al.1 (120111 ) calculate an apparent error rate of 0.10 mag, a bootstrap error of 
0.15 mag, and a .632 estimator of 0.13 mag based on optical-NIR light curves; for a more 
flexible model where th e Ry distributi o n diff ers in different Ay intervals they obtain a .632 
estimator of 0.11 mag. iFolatelli et al.l (|2010l ) calibrate absolute magnitudes in the redder 
optical bands r and i with training-set errors of 0.13 and 0.15 mag rms scatter and inferred 
intrinsic dispersions of 0.10 and 0.12 mag respectively. 

We compare the results of our methodology using the g-band calibration of the z = 0, 
gri filter set, with SALT2 in wavelength regions that correspond roughly to Johnson- 
Cousins B, V, and R. The combined subset dispersions and for the same set of supernovae 
are calculated to allow direct comparison with the preceding analysis. The results are 
wrms = 0.132 ± 0.008 mag and 4- fold cross-validation dispersion of 0.132 mag; the SALT2 
uncertainties in distance are large and result in an intrinsic dispersion consistent with zero. 
The linear-fit performs somewhat worse in reducing the total dispersions relative to this 
second-generation light-curve analyzer, although the Gaussian process method compares 
favorably. 

Most of the leverage in absolute magnitude comes from the first four parameters, as 
shown in Figure [T3j This contrasts with MLCS2k2 and SALT2, which describe supernovae 
with one free light-curve-shape and one free color parameter. Unlike those light curve 
fitters, many parameters are readily available in our analysis to account for color variance 
due to physical mechanisms within the supernova explosion and different distributions and 
properties of absorbing dust. 

The intrinsic dispersion is dependent on the blueshift of the filter set, a reasonable 
result con sidering th a t flux es from different restframe wavele ngth regions are measured. 
Note that Ijha et al.l ((2003, MLCS2k2) and iGuv et al.l (|20oj SALT2) each provide one 
redshift-independent intrinsic dispersion. As an example, the z-band of the z = 0.5 filter 
system has a dispersion of 0.114 mag and the i-band of the z = 0.75 filter system is 
0.131 mag. Given that the griz filters are not logarithmically distributed it is difficult 
to make direct inferences as to the exact cause of the different dispersions: whether local 
wavelength regions or the full wavelength range is playing the more important role in 
determining absolute magnitudes. It would be of interest to find correlations between the 
PCA coefficients of the z = 0.5 analysis with the excess residuals of the z = 0.75 analysis; 
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this would mean that systematic errors that depend on redshift are incurred due to the 
limited resolution and wavelength coverage of a fixed filter system. No such correlations 
are apparent with the available statistics. 

In the Gaussian process analysis, we restrict ourselves to using as parameters the co- 
efficients of the leading principal components that account for 95% of the dispersion. The 
reduction of the dimensionality of the parameter space was motivated by computational 
tractability and anticipating that it is overkill to use > 100 parameters to deduce absolute 
magnitudes. Indeed, our results do not change appreciably by using the ~ 30 parameters 
that describe 99% of the dispersion. Figure [13] hints that even our restricted parameter 
set is more than needed for this dataset. This validates our non-trivial expectation that 
the largest variances in light-curve shape and color be correlated with absolute magnitude. 
Further compression of information may be possible: In our methodology the PCA also 
accounts for the covariance in the light-curve prediction. This may be avoided by account- 
ing for measurement nois e while perfo rming the PCA, which is possible using expectation 



maximization techniques (|Bailevll2012l ). Other projections such as from independent com- 
ponent analysis and kernel PCA may prove more efficient in encapsulating the information, 
which could then be used in a Gaussian process. 

For a given supernova, we calculate the inferred peak absolute magnitude for each 
of the observer bands. These magnitudes and their residuals from the truth are highly 
correlated and so their resulting inferred distances are not independent. The same would 
be true if we had calibrated absolute magnitude at phases other than B maximum. The 
methodology presented optimizes the determination of absolute magnitude in a single band 
and phase. 

In an actual cosmology analysis, a customized light-curve model can be made for each 
redshift or on a grid sufficiently fine to make interpolation errors unimportant. For all 
redshifts, the distance moduli (difference between observed peak magnitude and inferred 
peak absolute magnitude) are based on a common set of supernovae and can be directly 
compared, at least to the same extent being done with current light-curve fitters. 



4. Future Outlook: Spectrophotometric Data 

The SNfactory generates spectrophotometry for objects at low redshift while the 
JDEM/ISWG design for a space telescope dark energy experiment called for spectrophoto- 
metric time series of high-redshift supernovae. In combination such data would bring the 
potential of providing improved distance measurements with lower statistical and system- 
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atic uncertainty by preserving information lost in broadband photometry. 



20()d : bailey et alJbood : IWang et alhood : iFolev fc Kaseilhoill : Ichotard et al 



3 magnitude ( 
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2011 
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2012h and may provide additional information not contained in broadband light curves that 
further standardize absolute magnitudes. Using the approach of $2] with each dispersion 
element serving as a filter, the variance in resolved features can be correlated with absolute 
magnitude. 

Supernoyae h ave varying magnitude dispersion depending on wavelength and phase 
(jJames et al that can be modeled as a Gaussian process. Measurements w,(i+ 2 )t,(i+ 2 )A 

for observer phase and wavelength (1 + z)t and (1 + z)\ are drawn as 

m ((i+z)t,(i+z)\) ~ GP (M(t, X;Pm) + 2.5 log {ji(z, pn)), k m (t, A, t', A'; -pk m ) + n m (t, A, t', A'; p„, m )) 

(11) 

With spectrophotometric data, the cosmological parameters = {Qm-, ^-de, u>0i w a} that 
determine the distance modulus [i predicted by theory can be fit simultaneously with the 
supernova parameters. 

Such an analysis using SNfactory data alone is of interest. Although insensitive to 
the cosmological parameters, the hyperparameters ~Pk m m particular provide insight into 
the interplay between wavelength and phase with the standard candle nature of Type la 
supernovae. This information aids in the planning of future missions and hardware and 
quantifies the systematic uncertainties potentially induced when using degraded data (e.g. 
broad-band photometry) only. 



5. Conclusions 

The regressed values from supernova light curves provide a large number of parameters 
that may correlate with absolute magnitude. Gaussian processes provide a convenient and 
parametrizable method of performing interpolations between measurements of light-curve 
shape and color and for using this information to predict absolute magnitudes. 

For a proof of concept, we present an example application showing how spectropho- 
tometry of low-redshift supernovae can be used to determine distances from broad-band 
photometry of high-redshift supernovae, such as will be obtained by DES and LSST. Using 
SNfactory data, we find dispersions in distance determinations as high as 0.44 mag cor- 
rected down to 0.04 mag within the training set and to as low as 0.09 mag when applied to 
a validation set. Some of this dispersion is attributable to the data themselves. Our proce- 
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dure is highly competitive with SALT2 and MLCS2k2, which give training-set dispersions 
of 0.12 and 0.14 respectively. 

The spectrophotometry of the SNfactory provides an excellent dataset for exploring 
supernova studies. Traditionally, broad-band light curves of z ~ supernovae are used 
to train light-curve fitters. The spectrophotometric templates allow the calibration of the 
absolute magnitude to be optimized for any redshift. With spectroscopy, a light-curve 
model can be tailor made for every object from an imaging survey. For planning of future 
supernova-specific imaging surveys, our methodology and SNfactory data can be used to 
determine optimal filter sets that minimize intrinsic dispersion over the targeted redshift 
range. 

There are directions for future development of the method and exploration of the data. 

• There are further improvements possible in the implementation of the methodology. 
A kernel other than the square exponential may better describe supernova light curves 
and their relationship with absolute magnitude. We did not perform an exhaustive 
search for the absolute maximum likelihoods in the parameter space since the local 
maxima found in this analysis were sufficient to establish our proof of concept. 



When constraining the problem to broad-band photometry, the full wavelength range 
and resolution available is not exploited. Analysis of intrinsic supernova properties 
within SNfactory and simultaneous analysis of low- and high-redshift samples are 
planned in future work. 

There is a correlation between supernov a absolute magnitudes and host-galaxy prop- 



erties after light-curve shape correction s (Kelly et al.ll2010l : ISullivan et al.ll2010l : lLampeitl et al 



2010i : lGupta et alJl201ll ; IChildress et al.ll2012l ). The search for such correlations with 



ga 
,2i 



the absolute magnitude residuals from our method is ongoing. 

• The observed variance of SNfactory measurements, and indeed of all supernova data, 
cannot be entirely attributed to photon noise. Gaussian processes provide a natural 
way to model and fit for external sources of photometric uncertainty. 
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A. Review of Gaussian Processes 

In this appendix we introduce Gaussian process and provide expressions for the likeli- 
hood and Normal probability distribution f unction used in this article. A broader descrip- 



tion of Gaussian processes can be found in iRasmussen fc Williams! ([200a) ■ 
A Gaussian process is denoted as 

/(x)~GP(m(x),t(x,x')). (Al) 

The expectation of the stochastic function /(x) is the m(x) with a model covariance be- 
tween any two points cov (/(x p ), f(x q )) = /c(x p ,x g ). Both m and k may be parameterized 
and the parameters of the latter are referred to as hyperparameters since they describe 
function scatter rather than the function itself. For a set of input points X+, the values of 
the function are drawn from a Normal distribution 

f*~Af(m(X*),if(X*,X*)), (A2) 

where the covariance matrix K has elements filled with all input pairs of fe(x*i,x*j). 

Since the probability distribution function is described by a Gaussian, several statistics 
are expressed analytically after algebraic manipulation. The likelihood that y, a set of n 
measurements of / at inputs X, with measurement covariance V described by a Gaussian 
process is written as 

logp(y|X) = -i(y - m{X)) T {K + V)-\y - m{X)) - I log \K + V\ - \ log 2vr. (A3) 
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A set of measurements and function values is drawn from a Normal distribution 



y 



m(X) 



K(X,X) + V K(X,X; 
K(X*,X) K(X*,X* 



(A4) 



The conditional distribution of the function values is Gaussian with expected mean 

f* = m(X) + K(X in X) [K(X, X) + V]' 1 (y - m(X)) (A5) 

and covariance 

cov (f*) = K{X*, X,) - K(X ir ,X) [K(X, X) + V]' 1 K(X, X,). (A6) 
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Table 4. True and Gaussian-process-inferred peak absolute magnitudes from validation 
using the z = filter set with a constant additive offset. In italics are supernovae that 

have data before two days pre-maximum with light curves that have at least eight points. 
In bold are supernovae that have data before two-days pre-maximum and at least eight 

points in all light curves. 
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10 


-14 


.84 


± 


0. 


03 


-14 


.73 


± 





.06 


-14 


12 


± 





04 


-14 


02 


± 


0. 


11 


73 


-15. 


24 


± 





.03 


-15 


.09 


± 


0. 


09 


-14 


99 


± 


0. 


04 


-14 


.86 


± 





.07 


-14 


38 


± 





04 


-14 


■ 21 


± 


0. 


10 


74 


-14 


74 


± 





.03 


-14 


.79 


± 





06 


-14 


.56 


± 





03 


-14 


.61 


± 





.04 


-13 


96 


± 





04 


-14 


02 


± 


0. 


07 


76 


-15. 


13 


± 





.04 


-15 


.19 


± 


0. 


06 


-14 


.90 


± 


0. 


06 


-14 


.94 


± 





.07 


-14 


18 


± 





.08 


-14 


22 


± 


0. 


06 


79 


-14 


82 


± 





04 


-14 


.99 


± 





11 


-14 


.64 


± 





04 


-14 


82 


± 





07 


-13 


94 


± 





05 


-14 


14 


± 





10 


81 


-15 


15 


± 





02 


-15 


.07 


± 





09 


-14 


.96 


± 





03 


-14 


89 


± 





07 


-14 


32 


± 





03 


-14 


23 


± 





06 


82 


-15 


20 


± 





03 


-15 


.20 


± 





09 


-14 


.90 


± 





02 


-14 


90 


± 





09 


-14 


23 


± 





04 


-14 


20 


± 





10 


84 


-15. 


03 


± 





.03 


-14 


.85 


± 


0. 


06 


-14 


.80 


± 


0. 


03 


-14 


.63 


± 





.04 


-H 


21 


± 





.03 


-H 


04 


± 


0. 


05 


85 


-14 


97 


± 





02 


-15 


.06 


± 





04 


-14 


.66 


± 





03 


-14 


74 


± 





04 


-14 


03 


± 





03 


-14 


12 


± 





05 


88 


-15 


44 


± 





03 


-15 


.23 


± 





06 


-15 


.16 


± 





03 


-14 


96 


± 





07 


-14 


56 


± 





04 


-14 


38 


± 





08 


89 


-15 


40 


± 





03 


-15 


.31 


± 





06 


-15 


.07 


± 





03 


-15 


01 


± 





04 


-14 


37 


± 





04 


-14 


26 


± 





07 


90 


-15 


00 


± 





03 


-15 


.21 


± 





08 


-14 


.74 


± 





03 


-14 


95 


± 





07 


-14 


04 


± 





03 


-14 


21 


± 





08 


91 


-15 


01 


± 





04 


-15 


.18 


± 





07 


-14 


.75 


± 





04 


-14 


90 


± 





05 


-14 


08 


± 





05 


-14 


22 


± 





07 


94 


-15 


02 


± 





03 


-15 


.20 


± 





05 


-14 


.71 


± 





03 


-14 


94 


± 





05 


-14 


01 


± 





04 


-14 


23 


± 





04 


95 


-15 


32 


± 





03 


-15 


.34 


± 





04 


-15 


.03 


± 





03 


-15 


05 


± 





03 


-14 


30 


± 





03 


-14 


31 


± 





05 


96 


-14 


■ 91 


± 





.03 


-15 


.13 


± 





08 


-14 


.67 


± 





04 


-14 


.88 


± 





.06 


-13 


87 


± 





.04 


-H 


08 


± 


0. 


10 


97 


-14 


.75 


± 





.03 


-14 


.85 


± 





07 


-14 


.63 


± 


0. 


02 


-14 


.71 


± 





.05 


-13 


85 


± 





.04 


-13 


■ 97 


± 


0. 


06 


98 


-14 


59 


± 





03 


-14 


.94 


± 





07 


-14 


.55 


± 





03 


-14 


85 


± 





04 


-13 


99 


± 





03 


-14 


38 


± 





07 


99 


-IS. 


96 


± 





03 


-14 


.21 


± 


0. 


03 


-14 


.16 


± 


0. 


03 


-14 


.34 


± 





.03 


-13 


63 


± 





03 


-IS 


89 


± 





.04 


100 


-15 


01 


± 





03 


-15 


.09 


± 





05 


-14 


77 


± 





03 


-14 


84 


± 





04 


-14 


07 


± 





03 


-14 


12 


± 





06 


102 


-14 


79 


± 





03 


-14 


76 


± 





09 


-14 


69 


± 





03 


-14 


64 


± 





07 


-14 


08 


± 





04 


-14 


07 


± 





10 


103 


-14 


92 


± 





03 


-14 


93 


± 





09 


-14 


69 


± 





04 


-14 


71 


± 





05 


-14 


15 


± 





04 


-14 


13 


± 





05 
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Table 4 — Continued 



SN # 






9 














r 












i 








true 




inferred 






true 






inferred 






true 




inferred 


105 


-15 


.12 ± 


03 


-15 


.06 ± 


10 


-14 


91 ± 





03 


-14 


.84 ± 


10 


-14 


.23 ± 


04 


-14 


18 ± 0.12 


107 


-15 


15 ± 


03 


-14 


.88 ± 


05 


-14 


89 ± 





03 


-14 


.66 ± 


06 


-14 


.27 ± 


.04 


-14 


00 ± 0.07 


108 


-15 


32 ± 


05 


-15 


.26 ± 


08 


-15 


00 ± 





06 


-14 


.93 ± 


06 


-14 


.31 ± 


06 


-14 


26 ± 0.07 


109 


-15 


18 ± 


03 


-15 


.18 ± 


09 


-14 


89 ± 





03 


-14 


.93 ± 


08 


-14 


.29 ± 


04 


-14 


33 ± 0.11 


111 


-15 


08 ± 


OS 


-15 


.13 ± 


.11 


-H 


76 ± 


0. 


OS 


-14 


.85 ± 


.08 


-14 


04 ± 


04 


-14 


08 ±0.00 


112 


-14 


71 ± 


04 


-14 


.66 ± 


08 


-14 


61 ± 





04 


-14 


59 ± 


05 


-14 


.10 ± 


06 


-14 


09 ± 0.05 


113 


-H 


77 ± 


.03 


-n 


.87 ± 


.07 


-14 


65 ± 





03 


-14 


. 73 ± 


.06 


-13 


97 ± 


04 


-14 


07 ± 0.07 


114 


-14 


86 ± 


03 


-14 


.95 ± 


09 


-14 


59 ± 





04 


-14 


.70 ± 


07 


-13 


92 ± 


05 


-14 


10 ± 0.12 


115 


-15 


52 ± 


03 


-15 


.30 ± 


11 


-15 


21 ± 





04 


-15 


02 ± 


09 


-14 


.65 ± 


04 


-14 


38 ± 0.11 


116 


-U 


53 ± 


.05 


-H 


. OH ± 


.07 


-U 


33 ± 


0. 


05 




.71 ±0 


.04 


-13 


81 ± 


.06 


-14 


20 ± 0.00 


117 


-15 


26 ± 


04 


-15 


.21 ± 


07 


-14 


98 ± 





04 


-14 


93 ± 


07 


-14 


21 ± 


06 


-14 


12 ± 0.11 


US 


-15 


32 ± 


04 


-15 


.19 ± 


09 


-15 


04 ± 





05 


-14 


93 ± 


07 


-14 


38 ± 


05 


-14 


25 ± 0.12 



Table 5. True and inferred peak absolute magnitudes from validation using the z = 0.25 

filter set. 



SN # true inferred true inferred true inferred true inferred 






-15 


07 


± 





07 


-15 


07 


± 





09 


-15 


08 


± 





04 


-15 


11 


± 





07 


-14 


94 


± 





04 


-14 


95 


± 





05 


-14 


08 


± 





05 


-14 


09 


± 





06 


1 


-14 


88 


± 





05 


-14 


85 


± 





09 


-15 


01 


± 





02 


-14 


94 


± 





06 


-14 


92 


± 





03 


-14 


90 


± 





05 


-14 


05 


± 





03 


-14 


05 


± 





05 


3 


-15 


11 


± 





05 


-15 


16 


± 





08 


-15 


12 


± 





02 


-15 


17 


± 





03 


-14 


98 


± 





02 


-15 


05 


± 





05 


-14 


11 


± 





03 


-14 


18 


± 





07 


4 


-14 


82 


± 





07 


-14 


67 


± 





13 


-14 


96 


± 





03 


-14 


84 


± 





09 


-14 


85 


± 





03 


-14 


70 


± 





09 


-14 


07 


± 





04 


-13 


96 


± 





09 


7 


-15 


.05 


± 





.07 


-15 


13 


± 





.11 


-15 


.04 


± 


0. 


03 


-15 


.19 


± 





07 


-14 


.93 


± 





.02 


-15 


05 


± 





.06 


-14 


02 


± 





03 


-14 


.11 


± 





.07 


8 


-15 


.02 


± 





.08 


-14 


80 


± 





15 


-15 


.05 


± 


0. 


04 


-14 


.88 


± 





.06 


-H 


.91 


± 





.04 


-14 


69 


± 





.04 


-14 


14 


± 





04 


-13 


.95 


± 





.05 


9 


-15 


02 


± 





05 


-15 


10 


± 





08 


-15 


02 


± 





02 


-15 


10 


± 





07 


-14 


87 


± 





03 


-14 


96 


± 





05 


-14 


00 


± 





03 


-14 


07 


± 





04 


10 


-15 


18 


± 





07 


-15 


11 


± 





10 


-15 


19 


± 





03 


-15 


13 


± 





06 


-15 


03 


± 





03 


-14 


96 


± 





05 


-14 


21 


± 





04 


-14 


14 


± 





06 


11 


-15 


.05 


± 





.10 


-14 


97 


± 





15 


-15 


12 


± 


0. 


04 


-15 


.03 


± 





07 


-14 


.99 


± 





.03 


-14 


92 


± 





.05 


-14 


12 


± 





05 


-14 


.01 


± 





.07 


12 


-15 


05 


± 





07 


-15 


02 


± 





13 


-15 


10 


± 





03 


-15 


09 


± 





15 


-14 


97 


± 





02 


-14 


90 


± 





11 


-14 


13 


± 





03 


-14 


09 


± 





13 


13 


-14 


29 


± 





06 


-14 


55 


± 





.09 


-14 


53 


± 





04 


-14 


.77 


± 





.06 


-14 


43 


± 





.04 


-14 


70 


± 





05 


-13 


70 


± 





04 


-13 


.95 


± 





.05 


14 


-14 


62 


± 





03 


-14 


60 


± 





10 


-14 


80 


± 





02 


-14 


76 


± 





.08 


-14 


75 


± 





02 


-14 


69 


± 





07 


-14 


01 


± 





02 


-14 


00 


± 





08 


16 


-14 


97 


± 





07 


-14 


86 


± 





12 


-15 


01 


± 





04 


-14 


84 


± 





09 


-14 


92 


± 





03 


-14 


78 


± 





08 


-14 


10 


± 





04 


-14 


00 


± 





09 


17 


-15 


06 


± 





07 


-14 


86 


± 





09 


-15 


23 


± 





04 


-14 


96 


± 





05 


-15 


01 


± 





03 


-14 


90 


± 





04 


-14 


24 


± 





05 


-14 


10 


± 





05 


21 


-14 


24 


± 





06 


-14 


24 


± 





08 


-14 


55 


± 





02 


-14 


56 


± 





04 


-14 


53 


± 





03 


-14 


53 


± 





04 


-13 


85 


± 





04 


-13 


84 


± 





03 


22 


-15 


38 


± 





06 


-15 


20 


± 





.09 


-15 


.30 


± 





03 


-15 


.10 


± 





.06 


-15 


.20 


± 





.03 


-14 


97 


± 





.06 


-14 


34 


± 





04 


-14 


.15 


± 





.07 


24 


-14 


50 


± 





07 


-14 


56 


± 





11 


-14 


71 


± 





02 


-14 


77 


± 





06 


-14 


68 


± 





03 


-14 


74 


± 





06 


-13 


90 


± 





04 


-13 


95 


± 





07 


25 


-15 


.05 


± 





10 


-15. 


03 


± 





.13 


-15 


.12 


± 


0. 


04 


-15 


10 


± 





.08 


-14 


.99 


± 





.05 


-14 


99 


± 





.07 


-14 


11 


± 





07 


-14 


.09 


± 





.08 


27 


-15 


37 


± 





07 


-15 


45 


± 





10 


-15 


28 


± 





04 


-15 


34 


± 





05 


-15 


09 


± 





.04 


-15 


16 


± 





03 


-14 


15 


± 





06 


-14 


23 


± 





05 


30 


-14 


90 


± 





06 


-14 


91 


± 





11 


-14 


99 


± 





02 


-14 


97 


± 





07 


-14 


88 


± 





03 


-14 


88 


± 





05 


-14 


11 


± 





04 


-14 


05 


± 





06 


33 


-15 


00 


± 





.08 


-15 


13 


± 





.12 


-15 


.03 


± 


0. 


03 


-15 


.13 


± 





.08 


-14 


.91 


± 





OS 


-15 


05 


± 





.06 


-14 


07 


± 





OS 


-14 


.19 


± 





.07 


36 


-15 


.21 


± 





.07 


-14 


80 


± 





.12 


-15 


.19 


± 


0. 


03 


-14 


.75 


± 





.05 


-15 


.15 


± 





.04 


-14 


75 


± 





.04 


-14 


62 


± 





04 


-14 


.06 


± 





04 


38 


-15 


29 


± 





06 


-15 


26 


± 





08 


-15 


25 


± 





03 


-15 


22 


± 





05 


-15 


07 


± 





03 


-15 


05 


± 





05 


-14 


22 


± 





04 


-14 


19 


± 





06 


39 


-15 


21 


± 





06 


-15 


26 


± 





10 


-15 


18 


± 





03 


-15 


23 


± 





05 


-14 


97 


± 





02 


-15 


09 


± 





04 


-14 


17 


± 





04 


-14 


25 


± 





05 


43 


-15 


.25 


± 





.08 


-15 


26 


± 





.13 


-15 


20 


± 





04 


-15 


.25 


± 





.06 


-15 


.07 


± 





.05 


-15 


10 


± 





06 


-14 


26 


± 





06 


-14 


.30 


± 





06 


45 


-15 


.37 


± 





.11 


-15 


30 


± 





16 


-15 


.32 


± 





06 


-15 


24 


± 





07 


-15 


.17 


± 





.05 


-15 


07 


± 





.04 


-14 


25 


± 





07 


-14 


.17 


± 





.07 


47 


-14 


.70 


± 





.07 


-14 


66 


± 





.08 


-14 


94 


± 


0. 


03 


-14 


.91 


± 





.05 


-14 


.87 


± 





.03 


-14 


84 


± 





05 


-14 


07 


± 





03 


-14 


.04 


± 





.04 


49 


-14 


.82 


± 





.05 


-14 


87 


± 





13 


-14 


94 


± 


0. 


03 


-14 


.99 


± 





.07 


-14 


.78 


± 





.04 


-14 


87 


± 





05 


-14 


02 


± 





04 


-14 


.09 


± 





.08 


54 


-15 


.20 


± 





.07 


-15 


09 


± 





.11 


-15 


.24 


± 


0. 


03 


-15 


.15 


± 





.06 


-15 


12 


± 





04 


-14 


99 


± 





.06 


-14 


14 


± 





06 


-H 


.02 


± 





.07 


55 


-14 


83 


± 





06 


-14 


96 


± 





10 


-14 


91 


± 





03 


-15 


02 


± 





10 


-14 


79 


± 





02 


-14 


92 


± 





08 


-13 


99 


± 





04 


-14 


11 


± 





08 


57 


-15 


26 


± 





04 


-15 


05 


± 





06 


-15 


29 


± 





02 


-15 


07 


± 





03 


-15 


12 


± 





02 


-14 


95 


± 





03 


-14 


30 


± 





03 


-14 


06 


± 





06 


59 


-14 


96 


± 





06 


-14 


95 


± 





12 


-14 


98 


± 





03 


-14 


98 


± 





07 


-14 


90 


± 





03 


-14 


90 


± 





07 


-14 


11 


± 





03 


-14 


09 


± 





07 


62 


-15 


12 


± 





06 


-15 


08 


± 





08 


-15 


10 


± 





03 


-15 


09 


± 





05 


-14 


94 


± 





03 


-14 


96 


± 





04 


-14 


19 


± 





04 


-14 


17 


± 





05 


63 


-14 


69 


± 





05 


-14 


63 


± 





09 


-14 


84 


± 





03 


-14 


77 


± 





06 


-14 


77 


± 





03 


-14 


73 


± 





06 


-14 


06 


± 





04 


-13 


97 


± 





06 


66 


-14 


60 


± 





07 


-14 


60 


± 





10 


-14 


77 


± 





03 


-14 


77 


± 





06 


-14 


72 


± 





03 


-14 


71 


± 





05 


-13 


90 


± 





04 


-13 


91 


± 





04 


67 


-15 


40 


± 





08 


-15 


25 


± 





12 


-15 


42 


± 





04 


-15 


19 


± 





.07 


-15 


26 


± 





04 


-15 


06 


± 





05 


-14 


43 


± 





06 


-14 


21 


± 





09 


69 


-14 


86 


± 





09 


-14 


86 


± 





13 


-14 


98 


± 





03 


-14 


93 


± 





07 


-14 


88 


± 





05 


-14 


86 


± 





07 


-14 


04 


± 





04 


-14 


05 


± 





06 


70 


-14 


63 


± 





07 


-14 


76 


± 





08 


-14 


73 


± 





04 


-14 


90 


± 





.07 


-14 


63 


± 





03 


-14 


75 


± 





07 


-13 


76 


± 





05 


-13 


95 


± 





05 


71 


-14 


.70 


± 





.06 


-14 


62 


± 





07 


-H 


.92 


± 


0. 


OS 


-14 


86 


± 





.05 


-14 


87 


± 





.03 


-14 


80 


± 





04 


-H 


05 


± 





04 


-13 


.98 


± 





.06 


73 


-15 


.03 


± 





.06 


-14 


90 


± 





H 


-15 


.13 


± 


0. 


03 


-14 


.93 


± 





.08 


-15 


.01 


± 





.03 


-14 


86 


± 





05 


-H 


18 


± 





04 


-14 


.00 


± 





.09 


74 


-14 


.48 


± 





07 


-14 


58 


± 





.12 


-14 


.63 


± 


0. 


02 


-H 


.75 


± 





.05 


-14 


.59 


± 





.03 


-H 


69 


± 





.05 


-13 


85 


± 





03 


-13 


.98 


± 





.05 


76 


-14 


.94 


± 





09 


-14 


99 


± 





10 


-15 


.01 


± 


0. 


06 


-15 


.07 


± 





.08 


-14 


92 


± 





.05 


-14 


98 


± 





.06 


-H 


09 


± 





08 


-14 


13 


± 





.08 


79 


-14 


49 


± 





07 


-14 


59 


± 





09 


-14 


71 


± 





04 


-14 


81 


± 





06 


-14 


68 


± 





04 


-14 


76 


± 





05 


-13 


88 


± 





04 


-13 


97 


± 





07 


81 


-14 


88 


± 





.04 


-14 


71 


± 





08 


-15 


06 


± 





02 


-14 


88 


± 





06 


-15 


00 


± 





02 


-14 


78 


± 





06 


-14 


16 


± 





03 


-14 


00 


± 





07 


82 


-15 


00 


± 





05 


-15 


04 


± 





08 


-15 


07 


± 





03 


-15 


13 


± 





05 


-14 


92 


± 





.04 


-15 


00 


± 





04 


-14 


11 


± 





03 


-14 


15 


± 





06 


84 


-H 


.87 


± 





.06 


-14 


83 


± 





.10 


-14 


.91 


± 


0. 


03 


-14 


.89 


± 





.05 


-14 


.82 







.03 


-14 


71 


± 





06 


-14 


07 


± 





03 


-14 


.06 


± 





04 


85 


-14 


80 


± 





05 


-14 


78 


± 





10 


-14 


83 


± 





02 


-14 


87 


± 





.04 


-14 


68 


± 





02 


-14 


65 


± 





03 


-13 


90 


± 





04 


-13 


91 


± 





03 


88 


-15 


37 


± 





05 


-15 


30 


± 





09 


-15 


34 


± 





03 


-15 


18 


± 





.07 


-15 


19 


± 





03 


-15 


06 


± 





06 


-14 


35 


± 





04 


-14 


26 


± 





09 


89 


-15 


25 


± 





07 


-15 


18 


± 





10 


-15 


25 


± 





03 


-15 


19 


± 





06 


-15 


11 


± 





03 


-15 


04 


± 





04 


-14 


25 


± 





04 


-14 


18 


± 





05 


90 


-14 


76 


± 





.04 


-14 


87 


± 





07 


-14 


87 


± 





03 


-14 


99 


± 





.04 


-14 


76 


± 





03 


-14 


90 


± 





05 


-13 


93 


± 





03 


-14 


09 


± 





05 


91 


-14 


79 


± 





06 


-15 


00 


± 





14 


-14 


87 


± 





03 


-15 


03 


± 





07 


-14 


76 


± 





04 


-14 


94 


± 





07 


-13 


98 


± 





05 


-14 


17 


± 





10 


94 


-14 


83 


± 





05 


-15 


05 


± 





09 


-14 


87 


± 





03 


-15 


08 


± 





.07 


-14 


75 


± 





03 


-14 


97 


± 





07 


-13 


92 


± 





03 


-14 


15 


± 





08 


95 


-15 


12 


± 





06 


-15 


13 


± 





08 


-15 


20 


± 





03 


-15 


22 


± 





04 


-15 


06 


± 





02 


-15 


07 


± 





04 


-14 


18 


± 





03 


-14 


20 


± 





06 


96 


-H 


.61 


± 





.07 


-H 


78 


± 





14 


-H 


.79 


± 


0. 


02 


-15 


.03 


± 





.09 


-14 


.70 


± 





.03 


-14 


87 


± 





11 


-13 


84 


± 





03 


-14 


.03 


± 





.12 


97 


-14 


.48 


± 





.07 


-14 


49 


± 





.11 


-14 


.65 


± 


0. 


02 


-14 


.68 


± 





.06 


-14 


.65 


± 





.02 


-H 


67 


± 





.07 


-13 


85 


± 





04 


-13 


.85 


± 





.09 


98 


-14 


.21 


± 





.05 


-H 


65 


± 





.07 


-14 


.56 


± 





02 


-14 


.98 


± 





.03 


-14 


.58 


± 





.03 


-14 


93 


± 





.03 


-13 


80 


± 





04 


-14 


.27 


± 





04 


99 


-13 


.36 


± 





06 


-13 


49 


± 





.08 


-13 


.99 


± 





02 


-14 


.17 


± 





.03 


-14 


.18 


± 





.02 


-14 


35 


± 





.03 


-13 


46 


± 





03 


-13 


.61 


± 





.02 


100 


-14 


82 


± 





05 


-15 


04 


± 





07 


-H 


.91 


± 





02 


-15 


11 


± 





.1)4 


-14 


.80 


- 





.03 


-15 


02 


± 





.04 


-13 


93 


± 





03 


-14 


16 


± 





05 


102 


-14 


55 


± 





06 


-14 


61 


± 





07 


-14 


71 


± 





03 


-14 


75 


± 





05 


-14 


71 


± 





03 


-14 


75 


± 





06 


-13 


90 


± 





04 


-13 


94 


± 





05 


103 


-14 


72 


± 





05 


-14 


75 


± 





08 


-14 


80 


± 





02 


-14 


89 


± 





06 


-14 


70 


± 





03 


-14 


76 


± 





05 


-13 


99 


± 





03 


-14 


03 


± 





05 
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Table 5 — Continued 



SN # 








9 














r 














i 














z 








true 






inferred 






true 






inferred 






true 






inferred 






true 






inferred 


105 


-14 


93 ± 





06 


-14 


.88 ± 


11 


-15 


.03 ± 





03 


-14 


99 ± 


07 


-14 


94 ± 





03 


-14 


82 ± 


05 


-14 


10± 





04 


-14 


.03 ± 0.08 


107 


-14 


98 ± 





06 


-14 


.94 ± 


08 


-15 


03 ± 





02 


-14 


87 ± 


06 


-14 


92 ± 





03 


-14 


84 ± 


06 


-14 


.15 ± 





03 


-14 


11 ± 0.05 


108 


-IS 


19 ± 





09 


-IS 


.11 ± 


16 


-IS 


.18 ± 


0. 


05 


-15. 


11 ±0 


.09 


-IS 


03 ± 


0. 


.05 


-14 


98 ± 


.09 


-14 


.20 ± 


0. 


07 


-14 


12 ± 0.08 


109 


-15 


00 ± 





06 


-15 


.05 ± 


11 


-15 


.05 ± 





02 


-15 


11 ± 


.08 


-14 


92 ± 





02 


-15 


01 ± 


07 


-14 


12 ± 





04 


-14 


.20 ± 0.08 


111 


-14 


93 ± 





06 


-IS 


25 ± 


09 




.94 ± 


0. 


03 


-15 


26 ±0 


.05 


-J4 


77 ± 





03 


-15 


13 ± 


.04 


-U 


. 00 ± 


0. 


04 




36 ± 0.05 


112 


-14 


36 ± 





08 


-14 


29 ± 


14 


-14 


63 ± 





04 


-14 


54 ± 


08 


-14 


63 ± 





03 


-14 


54 ± 


08 


-13 


93 ± 





04 


-13 


.87 ± 0.07 


113 


-U 


50 ± 





.07 




68 ± 


JJ 


-H 


72 ± 


0. 


03 


-14 


89 ± 


.07 


-J4 


68 ± 





03 


-J4 


86 ± 


05 


-13 


.84 ± 


0. 


04 




02 ± 0.05 


114 


-14 


63 ± 





.07 




77 ± 


13 


-H 


. 7,9 ± 


0. 


04 


-H 


91 ± 


.OS 


-14 


62 ± 


0. 


04 


-J4 


78 ±0 


.06 


-13 


.82 ± 


0. 


OS 


-13 


98 ±0.10 


115 


-15 


38 ± 





OS 


-15 


10 ± 


14 


-15 


.41 ± 





03 


-15 


11 ± 


09 


-15 


24 ± 





04 


-14 


97 ± 


09 


-14 


49 ± 





04 


-14 


18 ± 0.07 


116 


-14 


18 ± 





.08 


-H 


SO ± 


15 




.42 ± 


0. 


0^ 


-14 


77 ± 


.07 




35 ± 





05 


-14 


67 ± 


05 


-13 


.69 ± 


0. 


07 




00 ± 0.07 


117 


-15 


06 ± 





09 


-15 


.06 ± 


13 


-15 


.12 ± 





04 


-15 


14 ± 


07 


-15 


00 ± 





05 


-14 


99 ± 


04 


-14 


13 ± 





06 


-14 


13 ± 0.06 


118 


-IS 


09 ± 





.07 




.87 ± 


12 


-IS 


.18 ± 


0. 


04 


-14 


96 ± 


.05 


-15 


05 ± 





05 


-J4 


83 ± 


07 


-14 


.27 ± 


0. 


05 


-J4 


02 ± 0.05 
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Table 6. True and inferred peak absolute magnitudes from validation using the z = 0.50 

filter set. 



r i z 



SN # true inferred true inferred true inferred 






-15 


27 


± 





04 


-15 


.21 


± 





.08 


-14 


.98 


± 





.04 


-14 


90 


± 





06 


-14 


75 


± 





03 


-14 


69 







.05 


1 


-15 


14 


± 





03 


-15 


.22 


± 





07 


-14 


.91 


± 





.02 


-14 


83 


± 





06 


-14 


73 


± 





03 


-14 


69 


± 





.06 


3 


-15 


27 


± 





03 


-15 


27 


± 





07 


-15 


.01 


± 





03 


-15 


01 


± 





06 


-14 


79 


± 





03 


-14 


79 


± 





.06 


4 


-15 


11 


± 





04 


-14 


95 


± 





12 


-14 


.85 


± 





03 


-14 


72 


± 





09 


-14 


63 


± 





03 


-14 


50 


± 





.09 


6 


-15 


32 


± 





04 


-15 


.23 


± 





09 


-15 


03 


± 





04 


-14 


91 


± 





08 


-14 


82 


± 





03 


-14 


70 


± 





.07 


7 


-15 


22 


± 





.04 


-15 


.29 


± 





.08 


-U 


.96 


± 





.03 


-15 


.05 


± 





.06 


-14 


.75 


± 





.02 


-14 


.81 


± 





.07 


8 


-15 


23 


± 





05 


-15 


11 


± 





.08 


-14 


94 


± 





04 


-14 


81 


± 





06 


-14 


72 


± 





04 


-14 


58 


± 





.04 


9 


-15 


21 


± 





02 


-15 


.35 


± 





08 


-14 


92 


± 





03 


-15 


05 


± 





08 


-14 


68 


± 





02 


-14 


82 


± 





07 


10 


-15 


.27 


± 





03 


-15 


.10 


± 





.09 


-15 


10 


± 





.03 


-14 


.96 


± 





.08 


-14 


84 


± 





03 


-14 


76 


± 





.06 


12 


-15 


20 


± 





04 


-15 


15 


± 





11 


-15 


01 


± 





02 


-14 


96 


± 





10 


-14 


SO 


± 





02 


-14 


72 


± 





12 


13 


-14 


.60 


± 





03 


-14 


.77 


± 





.11 


-14 


42 


± 





.03 


-14 


61 


± 





07 


-14 


21 


± 





03 


-14 


46 


± 





.08 


14 


-14 


.89 


± 





.03 


-14 


.84 


± 





.11 


-14 


.70 


± 





.02 


-14 


.64 


± 





.11 


-14 


.53 


± 





.03 


-14 


41 


± 





.09 


16 


-15 


17 


± 





04 


-15 


14 


± 





13 


-14 


.92 


± 





03 


-14 


87 


± 





09 


-14 


70 


± 





03 


-14 


69 


± 





09 


17 


-15 


31 


± 





04 


-14 


.99 


± 





07 


-15 


.14 


± 





04 


-14 


76 


± 





06 


-14 


84 


± 





05 


-14 


55 


± 





.05 


19 


-19 


■ 91 


± 





.04 


-14 


.36 


± 





.07 


-13 


91 


± 





.03 


-14 


.37 


± 





.06 


-13. 


84 


± 





.04 


-14 


32 


± 





.05 


21 


-14 


55 


± 





03 


-14 


41 


± 





03 


-14 


44 


± 





02 


-14 


35 


± 





03 


-14 


31 


± 





03 


-14 


23 


± 





.03 


22 


-15 


49 


± 





03 


-15 


.29 


± 





.09 


-15 


.19 


± 





.02 


-15 


.01 


± 





.09 


-15 


00 


± 





03 


-14 


■ 79 


± 





.11 


24 


-14 


73 


± 





03 


-14 


.79 


± 





06 


-14 


.63 


± 





03 


-14 


66 


± 





05 


-14 


48 


± 





03 


-14 


50 


± 





.04 


27 


-15 


42 


± 





05 


-15 


.35 


± 





07 


-15 


.19 


± 





.04 


-15 


13 


± 





06 


-14 


91 


± 





04 


-14 


84 


± 





.04 


30 


-15 


07 


± 





03 


-15 


.27 


± 





09 


-14 


90 


± 





03 


-15 


08 


± 





10 


-14 


68 


± 





03 


-14 


87 


± 





08 


33 


-15 


15 


± 





.03 


-15 


.29 


± 





.08 


-U 


94 


± 





.03 


-15 


.07 


± 





.08 


-14 


71 


± 





03 


-14 


.83 


± 





.07 


36 


-15 


30 


± 





.04 


-15 


.02 


± 





.08 


-15 


.10 


± 





.03 


-14 


.82 


± 





.08 


-14 


95 


± 





.04 


-14 


67 


± 





.06 


38 


-15 


44 


± 





03 


-15 


40 


± 





06 


-15 


.15 


± 





03 


-15 


11 


± 





05 


-14 


88 


± 





03 


-14 


84 


± 





.06 


39 


-15 


36 


± 





03 


-15 


.50 


± 





07 


-15 


06 


± 





03 


-15 


18 


± 





06 


-14 


77 


± 





03 


-14 


87 


± 





.05 


45 


-15 


49 


± 





07 


-15 


.37 


± 





.11 


-15 


21 


± 





.05 


-15 


09 


± 





.06 


-14 


98 


± 





06 


-14 


87 


± 





.07 


47 


-14 


.94 


± 





.04 


-14 


.93 


± 





09 


-14 


.85 


± 





.03 


-14 


.86 


± 





.09 


-14 


67 


± 





.03 


-14 


66 


± 





.06 


49 


-15 


04 


± 





.03 


-15 


.15 


± 





.09 


-14 


.84 


± 





03 


-14 


.96 


± 





.07 


-14 


.59 


± 





.03 


-14 


.73 


± 





.05 


54 


-15. 


38 


± 





.05 


-15 


■ 24 


± 





.08 


-15 


.15 


± 





.04 


-15 


.01 


± 





.09 


-14 


.93 


± 





.04 


-14 


.78 


± 





.06 


55 


-15 


05 


± 





03 


-15 


.20 


± 





08 


-14 


.81 


± 





03 


-14 


91 


± 





OS 


-14 


60 


± 





03 


-14 


72 


± 





.07 


57 


-15 


40 


± 





03 


-15 


.23 


± 





04 


-15 


20 


± 





03 


-15 


05 


± 





03 


-14 


94 


± 





03 


-14 


80 


± 





.03 


59 


-15 


14 


± 





04 


-15 


.20 


± 





08 


-14 


.87 


± 





03 


-14 


96 


± 





08 


-14 


71 


± 





03 


-14 


73 


± 





.07 


62 


-15 


26 


± 





04 


-15 


.31 


± 





06 


-15 


.00 


± 





04 


-15 


06 


± 





05 


-14 


75 


± 





03 


-14 


81 


± 





.05 


63 


-14 


99 


± 





03 


-15 


.00 


± 





10 


-14 


.74 


± 





.02 


-14 


75 


± 





09 


-14 


56 


± 





03 


-14 


58 


± 





.10 


66 


-14 


88 


± 





04 


-14 


.94 


± 





.08 


-14 


.67 


± 





02 


-14 


74 


± 





08 


-14 


51 


± 





03 


-14 


55 


± 





.07 


67 


-15 


50 


± 





05 


-15 


.39 


± 





10 


-15 


.32 


± 





04 


-15 


19 


± 





11 


-15 


08 


± 





04 


-14 


97 


± 





.11 


69 


-15 


09 


± 





05 


-15 


.13 


± 





08 


-14 


.89 


± 





04 


-14 


93 


± 





07 


-14 


66 


± 





04 


-14 


72 


± 





.07 


70 


-14 


85 


± 





04 


-14 


.95 


± 





08 


-14 


.62 


± 





03 


-14 


74 


± 





07 


-14 


44 


± 





04 


-14 


51 


± 





.08 


71 


-15 


01 


± 





03 


-14 


.93 


± 





08 


-14 


.82 


± 





03 


-14 


78 


± 





06 


-14 


67 


± 





03 


-14 


62 


± 





.06 


73 


-15 


.23 


± 





.04 


-14 


.98 


± 





.10 


-15 


.04 


± 





03 


-U 


.79 


± 





.09 


-14 


82 


± 





.03 


-14 


.61 


± 





.08 


74 


-14 


75 


± 





.03 


-14 


.81 


± 





.10 


-14 


.51 


± 





03 


-14 


.62 


± 





.06 


-14 


36 


± 





.03 


-14 


43 


± 





.08 


76 


-15 


15 


± 





05 


-15 


.24 


± 





08 


-14 


95 


± 





.05 


-15 


05 


± 





07 


-14 


72 


± 





05 


-14 


84 


± 





06 


79 


-14 


81 


± 





04 


-14 


.90 


± 





09 


-14 


.62 


± 





03 


-14 


71 


± 





08 


-14 


47 


± 





04 


-14 


54 


± 





.07 


81 


-15 


10 


± 





02 


-14 


.90 


± 





07 


-14 


.98 


± 





02 


-14 


82 


± 





08 


-14 


81 


± 





03 


-14 


64 


± 





.06 


82 


-15 


20 


± 





03 


-15 


.27 


± 





.08 


-14 


.97 


± 





03 


-15 


01 


± 





07 


-14 


72 


± 





02 


-14 


78 


± 





.07 


84 


-15 


04 


± 





03 


-14 


.97 


± 





09 


-U 


.81 


± 





.02 


-14 


.75 


± 





.09 


-14 


61 


± 





03 


-14 


48 


± 





.06 


85 


-14 


99 


± 





03 


-15 


12 


± 





.05 


-14 


.72 


± 





03 


-14 


85 


± 





03 


-14 


48 


± 





03 


-14 


60 


± 





.03 


88 


-15 


42 


± 





03 


-15 


.30 


± 





.07 


-15 


.25 


± 





03 


-15 


15 


± 





08 


-15 


01 


± 





03 


-14 


91 


± 





.08 


89 


-15 


44 


± 





03 


-15 


.33 


± 





.08 


-15 


.15 


± 





.02 


-15 


02 


± 





07 


-14 


93 


± 





03 


-14 


79 


± 





.05 


90 


-15 


01 


± 





03 


-15 


.12 


± 





09 


-14 


.76 


± 





02 


-14 


86 


± 





07 


-14 


56 


± 





03 


-14 


69 


± 





.06 


91 


-15 


01 


± 





04 


-15 


.14 


± 





06 


-14 


.78 


± 





03 


-14 


92 


± 





06 


-14 


57 


± 





03 


-14 


70 


± 





.06 


94 


-15 


06 


± 





04 


-15 


22 


± 





.08 


-14 


77 


± 





.04 


-14 


92 


± 





08 


-14 


56 


± 





03 


-14 


71 


± 





.06 


95 


-15 


31 


± 





04 


-15 


.28 


± 





.05 


-15 


.10 


± 





03 


-15 


07 


± 





04 


-14 


88 


± 





03 


-14 


81 


± 





.04 


96 


-14 


.92 


± 





.05 


-15 


.16 


± 





.09 


-14 


.69 


± 





03 


-14 


.92 


± 





.07 


-14 


.50 


± 





.03 


-14 


.71 


± 





.06 


97 


-14 


79 


± 





03 


-14 


.86 


± 





09 


-14 


.53 


± 





.02 


-14 


63 


± 





07 


-14 


43 


± 





02 


-14 


49 


± 





.07 


98 


-14 


50 


± 





03 


-14 


76 


± 





.08 


-14 


46 


± 





.03 


-14 


.71 


± 





.04 


-14 


.38 


± 





.03 


-14 


M 


± 





.04 


99 


-IS. 


.78 


± 





.04 


-13 


.80 


± 





.06 


-13 


.92 


± 





.02 


-13 


.98 


± 





.04 


-13 


95 


± 





.03 


-13 


.98 


± 





.03 


100 


-15 


01 


± 





03 


-15 


17 


± 





05 


-14 


81 


± 





02 


-15 


00 


± 





04 


-14 


61 


± 





03 


-14 


78 


± 





.03 


102 


-14 


77 


± 





03 


-14 


79 


± 





06 


-14 


62 


± 





03 


-14 


59 


± 





05 


-14 


51 


± 





03 


-14 


48 


± 





06 


103 


-14 


93 


± 





04 


-14 


93 


± 





09 


-14 


68 


± 





03 


-14 


73 


± 





07 


-14 


49 


± 





03 


-14 


54 


± 





08 


105 


-15 


11 


± 





03 


-15 


07 


± 





09 


-14 


93 


± 





03 


-14 


91 


± 





07 


-14 


76 


± 





03 


-14 


73 


± 





.08 
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Table 6 — Continued 



SN # 






r 












i 










2 








true 




inferred 






true 




inferred 




true 




inferred 


107 


-15 


.19 ± 


04 


-15 


02 ± 


08 


-14 


91 ± 


03 


-14 


76 ± 


06 


-14.71 ± 


03 


-14 


.55 ± 0.07 


108 


-15 


37 ± 


.07 


-15. 


29 ± 


.11 


-i5. 


06 ± 


05 


-H 


.9.9 ± 


.07 


-14.85 ± 


.06 


-14 


76 ± 0.09 


109 


-15 


18 ± 


03 


-15 


21 ± 


08 


-14 


96 ± 


02 


-14 


99 ± 


09 


-14.73 ± 


03 


-14 


.74 ± 0.08 


111 


-IS. 


07 ± 


.04 


- 15 


19 ± 


10 




85 ± 


OS 




.93 ± 


.07 


-14.56 ± 


.03 




7^ ± 0.07 


112 


-14 


63 ± 


05 


-14 


53 ± 


10 


-14 


56 ± 


04 


-14 


.47 ± 


06 


-14.41 ± 


05 


-14 


36 ± 0.05 


113 


-H 


75 ± 


.05 




79 ± 


n 


-H 


6S ± 


OS 


-J4 


. 73 ± 


.09 


-14.49 ± 


.03 




57 ± 0.09 


114 


-H 


86 ± 


.04 


-J5. 


01 ± 


a 


-H 


62 ± 


OS 


-H 


. 74 ± 


.10 


± 


.04 


-14 


55 ± . 0* 


115 


-15 


53 ± 


.04 


-15 


26 ± 


ii 


-15 


31 ± 


03 


-15 


02 ± 


07 


-15.03 ± 


04 


-14 


.79 ± 0.07 


116 


-H 


± 


.05 


-H 


76' ± 


.08 


-H 


32 ±0 


05 


-J4 


. 61 ± 


.08 


-14.13 ± 


.05 


-i-i 


42 ± 0.06 


117 


-15 


.28 ± 


06 


-15 


20 ± 


08 


-15 


03 ± 


04 


-14 


96 ± 


05 


-14.82 ± 


04 


-14 


76 ± 0.06 


118 


-15 


33 ± 


05 


-15 


25 ± 


12 


-15 


07 ± 


04 


-14 


99 ± 


09 


-14.86 ± 


05 


-14 


76 ± 0.09 



Table 7. True and inferred peak absolute magnitudes from validation using the z = 

filter set. 







r 






i 






z 




SN # 


true 




inferred 


true 




inferred 


true 




inferred 






-14 


83 


± 





07 


-14 


83 


± 





10 


-15 


.18 


± 





04 


-15 


17 


± 





06 


-14 


73 


± 





04 


-14 


74 


± 





05 


1 


-14 


67 


± 





05 


-14 


58 


± 





09 


-15 


.04 


± 





02 


-14 


96 


± 





09 


-14 


68 


± 





02 


-14 


64 


± 





06 


3 


-14 


89 


± 





05 


-14 


89 


± 





10 


-15 


.18 


± 





03 


-15 


20 


± 





08 


-14 


78 


± 





03 


-14 


78 


± 





05 


4 


-14 


56 


± 





07 


-14 


34 


± 





13 


-15 


03 


± 





03 


-14 


81 


± 





11 


-14 


60 


± 





03 


-14 


40 


± 





10 


6 


-14 


81 


± 





07 


-14 


58 


± 





14 


-15 


22 


± 





03 


-14 


99 


± 





07 


-14 


79 


± 





03 


-14 


55 


± 





09 


7 


-14 


.89 


± 





.06 


-15 


.01 


± 





.12 


-15 


12 


± 





03 


-15 


.22 


± 





.09 


-14 


.70 


± 





03 


-14 


.82 


± 





.07 


8 


-14 


82 


± 





05 


-14 


65 


± 





07 


-15 


16 


± 





05 


-15 


01 


± 





10 


-14 


73 


± 





05 


-14 


55 


± 





05 


9 


-14 


82 


± 





05 


-14 


90 


± 





09 


-15 


09 


± 





02 


-15 


15 


± 





07 


-14 


69 


± 





02 


-14 


76 


± 





06 


10 


-15 


02 


± 





.06 


-14 


.94 


± 





11 


-15 


19 


± 





04 


-15 


.07 


± 





.06 


-14 


.87 


± 





03 


-14 


76 


± 





.06 


12 


-14 


86 


± 





06 


-14 


95 


± 





12 


-15 


.12 


± 





04 


-15 


18 


± 





10 


-14 


77 


± 





03 


-14 


87 


± 





10 


13 


-14 


02 


± 





.06 


-14 


.19 


± 





12 


-14 


54 


± 





04 


-14 


.73 


± 





.09 


-14 


18 


± 





04 


-14 


41 


± 





.09 


14 


-14 


33 


± 





05 


-14 


28 


± 





15 


-14 


.82 


± 





03 


-14 


80 


± 





16 


-14 


47 


± 





02 


-14 


38 


± 





12 


17 


-14 


86 


± 





07 


-14 


58 


± 





14 


-15 


.24 


± 





04 


-14 


95 


± 





11 


-14 


88 


± 





03 


-14 


58 


± 





12 


19 


-13. 


28 


± 





06 


-13 


.89 


± 





.11 


-13 


87 


± 





03 


-14 


49 


± 





05 


-IS 


68 


± 





03 


-14 


28 


± 





.05 


21 


-14 


02 


± 





04 


-14 


07 


± 





09 


-14 


.52 


± 





04 


-14 


58 


± 





07 


-14 


21 


± 





03 


-14 


31 


± 





06 


22 


-15 


■24 


± 





.07 


-15 


■ 15 


± 





■ 12 


-15 


.38 


± 





03 


-15 


.28 


± 





.09 


-14 


94 


± 





■ 03 


-14 


■ 86 


± 





■ 07 


24 


-14 


23 


± 





.06 


-14 


33 


± 


0. 


.10 


-14 


.66 


± 





04 


-14 


.77 


± 





09 


-14 


40 


± 





.03 


-14 


■ 45 


± 





.07 


27 


-15 


18 


± 





06 


-15 


18 


± 





09 


-15 


32 


± 





04 


-15 


29 


± 





05 


-14 


93 


± 





04 


-14 


89 


± 





04 


30 


-14 


74 


± 





06 


-14 


84 


± 





10 


-15 


.00 


± 





03 


-15 


15 


± 





07 


-14 


66 


± 





03 


-14 


82 


± 





09 


33 


-14 


.77 


± 





.05 


-14 


■ 91 


± 





08 


-15 


.05 


± 





OS 


-15 


18 


± 





.05 


-14 


70 


± 





.03 


-14 


81 


± 





.07 


36 


-15 


.06 


± 





.06 


-15 


.10 


± 





13 


-15 


.23 


± 





04 


-15 


.26 


± 





.08 


-14 


86 


± 





03 


-14 


.88 


± 





.08 


38 


-15 


10 


± 





06 


-15 


03 


± 





11 


-15 


35 


± 





03 


-15 


30 


± 





07 


-14 


91 


± 





03 


-14 


90 


± 





OS 


39 


-14 


.99 


± 





.05 


-14 


.98 


± 





10 


-15 


27 


± 





03 


-15 


.31 


± 





.07 


-14 


83 


± 





03 


-14 


.83 


± 





.05 


45 


-15. 


18 


± 





.08 


-15 


.10 


± 


0. 


11 


-15 


41 


± 


0. 


07 


-15 


.34 


± 





.08 


-14 


97 


± 





.06 


-14 


.90 


± 


0. 


04 


47 


-14 


44 


± 





.07 


-14 


.39 


± 


0. 


16 


-14 


89 


± 





03 


-14 


87 


± 





.10 


-14 


63 


± 





.03 


-14 


55 


± 


0. 


.09 


49 


-14 


.57 


± 





05 


-14 


47 


± 


0. 


14 


-15 


00 


± 





04 


-14 


90 


± 





.14 


-14 


61 


± 





03 


-14 


54 


± 


0. 


.12 


54 


-14 


96 


± 





07 


-14 


81 


± 





14 


-15 


.28 


± 





05 


-15 


13 


± 





.11 


-14 


92 


± 





04 


-14 


.75 


± 





.07 


55 


-14 


60 


± 





04 


-14 


63 


± 





11 


-14 


.98 


± 





03 


-14 


99 


± 





07 


-14 


57 


± 





03 


-14 


58 


± 





08 


57 


-15 


14 


± 





04 


-14 


96 


± 





11 


-15 


.31 


± 





03 


-15 


21 


± 





07 


-14 


96 


± 





03 


-14 


83 


± 





06 


59 


-14 


73 


± 





06 


-14 


76 


± 





09 


-15 


.04 


± 





03 


-15 


06 


± 





06 


-14 


63 


± 





03 


-14 


70 


± 





05 


62 


-14 


90 


± 





05 


-14 


88 


± 





09 


-15 


20 


± 





03 


-15 


21 


± 





08 


-14 


77 


± 





03 


-14 


76 


± 





07 


63 


-14 


■ 36 


± 





.05 


-14 


■24 


± 





10 


-H 


.88 


± 





04 


-H 


.73 


± 





.11 


-14 


■ 51 


± 





.03 


-14 


.38 


± 


0. 


.09 


66 


-14 


33 


± 





.06 


-14 


.35 


± 





14 


-14 


.78 


± 





04 


-14 


.82 


± 





.09 


-14 


45 


± 





.03 


-14 


42 


± 


0. 


.10 


67 


-15 


25 


± 





06 


-15 


14 


± 





10 


-15 


42 


± 





05 


-15 


30 


± 





09 


-15 


06 


± 





04 


-14 


94 


± 





06 


69 


-14 


58 


± 





05 


-14 


41 


± 





13 


-15 


.03 


± 





05 


-14 


90 


± 





11 


-14 


65 


± 





04 


-14 


53 


± 





12 


70 


-14 


.30 


± 





.07 


-14 


43 


± 





14 


-14 


.80 


± 





04 


-H 


.99 


± 





.12 


-14 


39 


± 





.03 


-14 


.55 


± 


0. 


.11 


71 


-14 


40 


± 





.05 


-14 


23 


± 





.08 


-14 


.95 


± 





OS 


-14 


74 


± 





.11 


-14 


.59 


± 





.03 


-14 


40 


± 


0. 


.09 


73 


-14 


84 


± 





.05 


-14 


56 


± 





11 


-15 


.13 


± 





04 


-14 


.83 


± 





.09 


-14 


80 


± 





.03 


-14 


■ 54 


± 


0. 


.07 


74 


-14 


23 


± 





.05 


-14 


■ 58 


± 





IS 


-14 


.68 


± 





03 


-14 


.99 


± 





.12 


-14 


28 


± 





.03 


-14 


64 


± 


0. 


.08 


76 


-14 


70 


± 





.06 


-14 


68 


± 





.08 


-15 


.04 


± 





04 


-15 


.02 


± 





.05 


-14 


.70 


± 





.05 


-14 


■ 66 


± 


0. 


.05 


79 


-H 


12 


± 





.07 


-14 


10 


± 





.11 


-U 


■ 74 


± 





OS 


-14 


.71 


± 





.06 


-H 


.39 


± 





.04 


-14 


.35 


± 


0. 


06 


81 


-14 


67 


± 





04 


-14 


49 


± 





09 


-15 


.03 


± 





03 


-14 


78 


± 





09 


-14 


75 


± 





02 


-14 


53 


± 





06 


82 


-14 


77 


± 





05 


-14 


79 


± 





08 


-15 


13 


± 





03 


-15 


13 


± 





09 


-14 


73 


± 





03 


-14 


73 


± 





06 


84 


-H 


.66 


± 





.04 


-14 


72 


± 





IS 


-H 


.97 


± 





OS 


-15 


.06 


± 





.11 


-H 


.58 


± 





.03 


-14 


.60 


± 


0. 


■ 08 


85 


-14 


58 


± 





05 


-14 


78 


± 





05 


-14 


.92 


± 





03 


-15 


14 


± 





03 


-14 


47 


± 





03 


-14 


68 


± 





04 


88 


-15 


27 


± 





06 


-15 


24 


± 





11 


-15 


36 


± 





03 


-15 


33 


± 





09 


-15 


00 


± 





03 


-15 


00 


± 





04 


89 


-15 


07 


± 





06 


-14 


98 


± 





10 


-15 


.34 


± 





03 


-15 


25 


± 





08 


-14 


91 


± 





03 


-14 


87 


± 





07 


90 


-14 


51 


± 





.07 


-14 


59 


± 





14 


-14 


.93 


± 





03 


-14 


99 


± 





10 


-14 


53 


± 





04 


-14 


59 


± 





09 


91 


-H 


.51 


± 





.07 


-14 


58 


± 





.09 


-H 


.94 


± 





04 


-15 


.02 


± 





.06 


-H 


.54 


± 





.04 


-14 


.61 


± 


0. 


04 


94 


-14 


53 


± 





05 


-14 


74 


± 





13 


-14 


.95 


± 





03 


-15 


20 


± 





10 


-14 


54 


± 





03 


-14 


74 


± 





09 


95 


-14 


91 


± 





05 


-14 


87 


± 





10 


-15 


.22 


± 





03 


-15 


18 


± 





07 


-14 


86 


± 





03 


-14 


82 


± 





05 


96 


-H 


33 


± 





.06 


-14 


.54 


± 





13 


-14 


.83 


± 


0. 


OS 


-15 


.04 


± 





.13 


-H 


47 


± 





.04 


-14 


.65 


± 


0. 


11 


97 


-14 


10 


± 





.05 


-14 


.16 


± 





13 


-14 


.68 


± 


0. 


OS 


-14 


.72 


± 





.13 


-14 


.32 


± 





.02 


-14 


41 


± 


0. 


■ 10 


98 


-14 


17 


± 





.04 


-14 


.79 


± 





.09 


-14 


44 


± 





.03 


-15 


.13 


± 





.06 


-14 


23 


± 





.02 


-14 


.76 


± 


0. 


05 


99 


-is. 


00 


± 





.06 


-13. 


43 


± 





.10 


-13 


.75 


± 





.03 


-14 


.20 


± 





.05 


-is 


.69 


± 





.03 


-14 


04 


± 


0. 


05 


100 


-14 


52 


± 





.07 


-H 


60 


± 


0. 


11 


-H 


.91 


± 


0. 


OS 


-15 


.01 


± 





.06 


-14 


57 


± 





.03 


-14 


68 


± 


0. 


06 


102 


-14 


34 


± 





04 


-14 


62 


± 





11 


-14 


69 


± 





03 


-14 


90 


± 





07 


-14 


38 


± 





03 


-14 


64 


± 





09 


103 


-14 


50 


± 





05 


-14 


55 


± 





12 


-14 


86 


± 





03 


-14 


91 


± 





07 


-14 


44 


± 





03 


-14 


48 


± 





09 


105 


-14 


72 


± 





05 


-14 


80 


± 





11 


-15 


03 


± 





03 


-15 


12 


± 





07 


-14 


68 


± 





03 


-14 


74 


± 





11 


107 


-14 


75 


± 





05 


-14 


67 


± 





10 


-15 


10 


± 





03 


-15 


01 


± 





07 


-14 


66 


± 





03 


-14 


59 


± 





07 


108 


-14 


92 


± 





.08 


-H 


83 


± 





H 


-15 


.29 


± 





04 


-15 


.21 


± 





.10 


-H 


.83 


± 





.06 


-14 


75 


± 





09 



-47- 



Table 7 — Continued 



SN # 


r 

true 


inferred 


i 

true 


inferred 


z 

true 


inferred 


109 


-14.85 ± 0.06 


-14.97 ± 0.13 


-15 .10 ± 0.04 


-15.20 ± 0.10 


-14.72 ± 0.03 


-14.82 ± 0.09 


111 


-14.54 ± 0.0,5 


-14.51 ± 0.10 


-14.98 ± 0.04 


± 0.07 


-14.62 ± 0.03 


-14. 60 ± 0.06 


112 


-I4.IO ± 0.07 


-14.17 ± 0. J,5 


-14.59 ± 0.05 


-14.72 ± O.JO 


-14.3S ± O.O^ 


- J^46 ± 0.0* 


113 


-14.21 ± 0.06 


-14.SS ± O.J J 


-14.64 ± 0.04 


-14.70 ± 0.0,9 


-14.41 ± 0.03 


- J^.SO ± 0.06 


114 


-14.S4 ± 0.07 


- J^.53 ± O.J 6 


-14. 81 ± 0.0-i 


-15.06 ± 0. J2 


-14.39 ± O.O^ 


-14.61 ± 0.0.9 


118 


-15.18 ± 0.07 


-14.88 ± 0.12 


-15.44 ± 0.04 


-15.15 ± 0.09 


-15.06 ± 0.03 


-14.82 ± 0.10 


116 


-1S.91 ± 0.06 


- J^.25 ± 0.12 


± 0.05 


-14.78 ± 0.0S 


-14.10 ± 0.0,5 


-14.45 ± 0.06 


117 


-14.88 ± 0.08 


-14.84 ± 0.14 


-15.17± 0.05 


-15.12 ± 0.09 


-14.79 ± 0.03 


-14.76 ± 0.08 


118 


-14.87 ± 0.06 


-14.57 ± 0.15 


-15.24 ± 0-05 


-14.96 ± 0.0,9 


- 14.83 ± 0.0,5 


-14.53 ± 0.08 



